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A nonlinear  theory  of  elastic  sandwich  plates  has 
been  developed  from  the  Modified  Hel 1 inger-Rei ssner  prin- 
ciple in  three-dimensional  elasticity  using  a variational 
method.  A consistent  method  of  reducing  in  steps  the  most 
general  three-dimensional  mixed  boundary  value  problem  to, 
first,  a two-dimensional  mixed  boundary  value  problem,  and 
then  to  a single  layer  or  multi-layer  plate  problem  by 
assuming  different  displacement  functions,  has  been  pre- 
sented. A complete  set  of  equilibrium  equations,  boundary 
conditions,  compatibility  equations  and  strain  functions 
in  terms  of  displacement  functions  is  derived  in  a unified 
approach  from  the  variational  principle  mentioned  above, 
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whereas  in  the  earlier  theories  of  other  authors,  the 
compatibility  equations  and  the  relations  between  strain 
functions  and  displacement  functions  had  to  be  considered 
separately  and  a priori. 

The  existing  analysis  of  an  elastic  sandwich  plate 
subjected  to  combined  axial  and  transverse  loads  and  the 
corresponding  buckling  problem  are  shown  to  be  special 
cases  of  the  general  analysis  presented  here. 

In  the  second  part  of  this  dissertation,  a known 
solution  of  the  problem  of  a rectangular  sandwich  plate, 
simply  supported  on  all  sides  and  subjected  to  arbitrary 
transverse  loads  on  both  sur faces  a n o uniformly  distributed 
uniaxial  loads  on  two  opposite  edges,  has  been  used  to  get 
the  solution  of  a similar  problem  of  a sandwich  plate  made 
of  viscoelastic  materials,  using  the  correspondence  prin- 
ciple of  viscoelasticity  theory.  Expressions  for  time 
dependent  antisymmetric  and  symmetric  deflections  and  time 
dependent  critical  values  of  the  compressive  axial  loads 
for  both  antisymmetric  and  symmetric  modes  of  buckling  have 
been  derived.  Calculations  have  been  done  and  results  have 
been  shown  for  some  specific  problems. 
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CHAPTER  1 


INTRODUCTION 


General  Background 


A sandwich  plate  is  the  plate  form  of  structural 
sandwich  construction  which  is  defined  [1]*  by  the  Ameri- 
can Society  for  Testing  and  Materials  as: 

A laminar  construction  comprised 
of  alternate  dissimilar  simple  or  com- 
posite materials  assembled  and  intimately 
fixed  in  relation  to  each  other  so  as  to 
use  the  properties  of  each  to  attain  spe- 
cific structural  advantages  for  the  whole 
assembly  (p.  1.3). 

A sandwich  structure  generally  consists  of  rela- 
tively thin  external  members  called  facings  separated  by 
and  bonded  to  a relatively  thick  low  density  internal  mem- 
ber called  a core.  Usually  the  facings  are  composed  of 
high  strength  materials  and  the  core  is  made  of  materials 
having  low  average  density  and  strength. 

Sandwich  structures  exhibit  the  following  advantages 
high  ratio  of  bending  stiffness  to  weight,  resistance  to 


*Numbers  in  brackets  refer  to  Bibliography. 
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fatigue,  sound  damping,  thermal  i nsul ati on  and  aerodynam- 
ically  smooth  surfaces. 

Continuous  improvements  in  the  techniques  and  tech- 
nology of  adhesive  bonding  have  made  the  adhesive-bonded 
sandwich  structures  practical  in  an  ever-increasing  range 

of  applications.  Examples  are: 

Aircraft:  Control  surfaces,  control  tabs,  bulk- 

heads, helicopter  blades,  flooring,  doors,  engine 
cowlings,  wing  panels,  fairings,  fan  air  ducts, 
etc  . 

Missile  an  d spacecraft : Fins  and  control  surfaces, 

tanks,  structural  shells,  heat  shields,  intertanK 
structures  , etc . 

El ectronics  and  communication  : Antenna  reflectors, 

radomes,  electronic  packaging,  etc. 

Building:  Pre-fabri cated  shelters,  curtain  walls, 

partitions,  doors,  etc. 

Miscellaneous:  Desk  top  scaffolding,  house  trailer 

flooring,  boats,  cargo  containers,  etc. 

Depending  on  the  conditions  under  which  it  is  used, 
the  material  of  the  faces  may  be  an  aluminum  alloy,  rein- 
forced plastic,  titanium,  heat-resistant  steel,  etc.  A 
very  popular  type  of  core  now  in  use  is  the  'honeycomb 
core,  which  consists  of  thin  foils  in  the  form  of  hexagonal 
cells  perpendicular  to  the  faces.  Other  types  of  core  are 
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corrugated  sheet,  expanded  materials  like  cellulose  ace- 
tate, synthetic  rubber,  balsa  wood,  etc.  The  adhesive 
material  used  for  bonding  the  layers  to  the  core  may  be 
Nitrile  Phenolic,  Vinyl  Phenolic,  Epoxy  Phenolic,  Unmodi- 
fied Epoxy,  Epoxy  Polymide,  etc.  The  dimensions  of  a sand- 
wich construction  vary  widely  from  one  case  to  another.  A 
typical  ratio  of  core  thickness  to  the  face  thickness  may 
be  20. 

Although  there  were  some  earlier  applications,  sand- 
wich construction  aroused  great  interest  only  after  the 
Second  World  War  when  1 ami na r- f 1 ow  surface  was  considered 
to  be  an  extremely  desirable  design  feature  for  high  speed 
aircraft  and  when  it  became  known  that  the  famous  Mosquito 
airplane  was  constructed  in  this  manner.  Therefore,  the 
sandwich  construction  is  a comparatively  new  branch  of 
structural  elements  with  a history  of  about  30  years.  But 
because  of  its  wide  range  of  applicability,  it  has  attrac- 
ted the  attention  of  hundreds  of  investigators  from  all 
over  the  world,  the  sum  total  of  whose  contributions  has 
established  its  technology  on  a highly  refined  and  ever- 
improving  theoretical  basis.  The  total  number  of  articles 
published  in  this  field  (well  over  1,000  articles  and  more 
than  10  books)  also  signifies  the  importance  of  sandwich 
constructi on . 

The  theory  and  analysis  of  the  sandwich  plate  are 
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much  more  involved  than  and  different  from  those  of  ordinary 
single  layer  plates  because  the  response  of  sandwich  con- 
struction to  various  inputs  is  determined  to  a large  ex- 
tent by  the  geometrical  and  material  properties  of  the  core 
relative  to  the  facings.  The  single  layer  plate  can  be  con- 
sidered to  be  a special  case  of  the  sandwich  plate  and  the 
analysis  of  the  former  can  be  derived  as  a special  case  from5 
and  be  used  as  a simple  check  on, the  analysis  of  the  latter. 

A general  survey  of  the  extensive  literature  on 
sandwich  plates  provides  the  history  and  the  present  state 
of  the  art  which  cover  the  effect  of  transverse  shear  de- 
formation on  the  small  as  well  as  large  deflections  of  sand- 
wich plates,  the  locai  instability  due  to  wrinkling  of  the 
facings,  the  intracell  buckling  cf  honeycomb  sandwiches, 
the  general  instability  of  sandwich  structures,  the  thermal 
effects  as  well  as  the  problems  of  dynamic  response  including 
the  effects  of  various  inertia  terms,  damping,  and  large 
amplitudes.  A few  reports  on  viscoelastic  behavior  and 
wave  propagation  problems  of  sandwich  plates  are  also  avail- 
able. 

In  the  following,  an  attempt  is  made  to  recount 
the  most  prominent  phases  of  development  of  the  theory  of 
sandwich  plates  while,  for  fuller  surveys,  it  is  suggested 
that  the  reader  see  the  comprehensive  bibliographies  [2], 

[3]  and  the  review  articles  [4]  and  [5]. 


Hi s t o r y and  Present  State  of  the  Art 


Sma 1 1 Deflection  Theory 

The  early  theory,  which  accounts  for  the  transverse 
shear  effect  in  the  core  by  assuming  that  a linear  element 
initially  straight  and  normal  to  the  middle  plane  of  the 
core  will  remain  straight  after  deformation  but  will  de- 
viate from  the  normal  to  the  deformed  middle  plane  by  an 


amount  represented  by  a parameter  (wh i ch  depends  on  the  elas 
tic  properties  of  the  material),  was  developed  by  a number 


of  authors,  e.g.  [6],  [7]. 

B a t d o r f [S]  andReissner[9] 


The  later  theory  of  Libove  and 

has  led  to  the  deflection  equa- 


tions of  sandwich,  plates  on  the  basis  of  a simplified  model. 
Reissner  [10]  noted  that,  for  such  simplifying  assumptions, 
the  equilibrium  equations  of  a sandwich  plate  have  the 
same  form,  except  for  constants,  as  his  equations  for  homo- 
geneous plates  with  transverse  shear  deformation  taken  into 
account.  Various  effects  previously  omitted  have  been  in- 
corporated into  the  analysis  of  Eringen  [11]  who  generalized 
the  variational  approach  used  by  Hoff  [12].  A comparatively 
recent  attempt  by  Mushtari  [13]  has  resulted  in  considerable 


refinement  of  earlier  theories. 

Thermal  effects  in  sandwich  plates  have  been  inves- 
tigated by  Egcioglu  [14],  among  others,  with  the  help  of 
suitable  transf ormati ons  of  the  variables  in  order  to  obtain 
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two  i ndependent  systems  of  differential  equations.  More 
recent  works  including  [15]  and  [16]  show  minor  improve- 
ments of  the  theory. 

Large  Deflections 

Reissner  [17]  employed  a simplified  model  consist- 
ing of  two  facings  acting  as  membranes  and  a core  resisting 
transverse  shear  and  normal  stresses  only  for  deriving  a 
set  of  equations  governing  the  large  deflection  of  sand- 
wich plates.  The  principle  of  complementary  energy  has  been 
used  by  Teichmann  and  Wang  [18]  for  curved  sandwich  plates 
and  cylinders  in  deriving  non-linear  stress  resultant-dis- 
placement relations. 

Several  contributions  to  the  large  deflection  analy- 
sis of  sandwich  structures  have  been  published  in  the  Sov- 
iet Union  over  the  past  fifteen  years.  Many  of  these  have 
been  summarized  by  Habip  [4 ] . They  contain  analyses  of 
orthotropic  cores  using  both  the  variational  and  the  clas- 
sical methods. 

Recently  Ebcioglu  [19],  [20]  has  given  a non-linear 
theory  of  sandwich  panels  in  the  reference  state  using 
Hamilton's  principle.  He  obtained  the  equations  of  motion 
and  the  boundary  conditions  for  twelve  displacement  func- 
tions. The  analysis  has  included  orthotropic  materials. 
Numerical  solutions  of  specific  problems  have  been  worked 
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out  by  many  authors,  e.g.  [21]. 

General  Buckling 

Of  all  the  earlier  work  done  in  this  field,  the 
one  by  Eringen  [22]  stands  as  the  most  prominent,  and  it 
gives  solutions  for  both  the  antisymmetric  and  symmetric 
modes  of  buckling.  The  work  of  Hoff  [12]  is  also  worth 
mentioning  in  this  context.  However,  all  these  works  were 
associated  with  simple  boundary  conditions.  Nardo  [23] 
obtained  an  exact  solution  for  a rectangular  sandwich  plate 
with  loaded  edges  clamped  and  the  other  two  edges  simply 
supported.  The  energy  method  has  been  used  by  Ericksen 
and  March  and  they  published  a series  of  papers  [24],  [25], 
[26]  on  the  solution  of  problems  for  a number  of  useful 
boundary  conditions.  In  their  solutions,  the  materials  of 
both  the  faces  and  the  core  may  be  orthotropic.  An  im- 
provement of  the  Ericksen-March  method  is  given  in  [27], 
which  gives  solution  for  the  buckling  by  uniaxial  compres- 
sive load.  Moskalanko  [28]  obtained  linearized  equations 
of  non  linear  elasticity  and  also  found  the  components  of 
the  stress  tensor  perturbations  across  the  plate. 

Among  the  recent  reports  on  experimental  determin- 
ation of  the  critical  stresses  for  buckling,  the  one  by 
Zender  [29]  gives  results  for  a large  number  of  cases.  He 
has  compared  his  results  with  the  theoretical  results  of 
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orthotropic  plates.  The  comparison  shows  that  the  open- 
face  truss-core  sandwich  panel  is  substantially  less  effi- 
cient than  indicated  by  theory.  The  inefficiency  was  per- 
haps the  result  of  reduced  stiffness  near  the  edges. 

Numerical  results  have  been  obtained  by  Benson  [30] 
arid  Landgren  and  Salama  [31]  on  the  solution  of  buckling 
problems.  Benson  considered  both  sandwich  plates  and  col- 
umns with  the  effects  of  shear  and  normal  stresses  in  the 
core  which  is  orthotropic.  Using  variational  methods,  he 
obtained  differential  equations  governing  the  phenomena  of 
both  general  buckling  and  face  wrinkling.  Landgren  and 
Salama  used  finite  element  methods  to  solve  the  problem  of 
buckling  of  multilayer  plates.  Ebcioglu  and  Kim  [32]  mod- 
ified some  of  the  assumptions  made  by  earlier  authors  and 
solved  both  the  antisymmetric  and  the  symmetric  problems 
of  buckling  of  rectangular  sandwich  plates  making  use  of 
the  continuity  relations  at  the  interfaces  between  layers. 

Local  Buckling 

Of  the  earlier  reports  on  the  problem  of  local 
buckling,  Norris  [33]  and  Eringen  [22]  presented  relatively 
simple  methods  for  finding  the  wrinkling  stresses  in  the 
facing  of  an  orthotropic  sandwich  construction.  Yusuff's 
[34]  theory  took  into  consideration  the  initial  imperfec- 
tions of  sandwich  plates.  A different  type  of  local 
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instability  due  to  the  particular  geometry  of  the  core, 
such  as  honeycomb,  is  intracell  buckling  of  the  facings. 
Weikel  and  Kobayashi  [35]  have  reported  some  results  ob- 
tained under  various  simplifying  assumptions  and  compared 
them  with  the  earlier  attempts  of  Norris  and  Kommers  [36]. 
Recently,  Benson  [30]  has  presented  both  theory  and  numer- 
ical results  based  on  variational  principles.  They  also 
considered  general  buckling  as  mentioned  earlier. 

Apart  from  the  areas  described  above,  vibration 
problem  of  sandwich  plates  has  been  treated  by  several 
authors  including  Yu  [37],  [38]  who  has  investigated  the  sim- 
ple-shear modes  of  vibration  of  sandwich  plates  on  the 
basis  of  a comprehensive  sandwich  plate  theory  and  also 
the  nonlinear  flexural  vibrations  of  sandwich  plates. 

The  high  frequency  extensional  vibrations  of  s a n d w i c h 
plates  have  been  investigated  by  Chu  [39],  The  wave  propa- 
gation problem  in  sandwich  plates  has  been  investigated  by 
Armenakas  [40],  among  others. 

Previous  Work  Related  to  and  the 
Objectives  of  This  Dissertation 

As  mentioned  earlier  in  this  section,  Ebcioglu 
[19],  [20]  obtained  a complete  set  of  equations  of  motion 
and  boundary  conditions  from  Hamilton's  principle  using 
the  variational  method.  However,  the  equations  relating 
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various  stress  resultants  to  strain  functions  and  the 
relations  between  strain  functions  and  displacement  func- 
tions had  to  be  considered  separately  andapriori.  Habip 
[41]  used  a different  procedure.  He  started  with  the  Modi- 
fied Hellinger-Reissner's  Principle  and  obtained  the  rela- 
tions between  stress  resultants  and  strain  functions,  and 
between  strain  functions  and  displacement  functions,  in 
addition  to  the  equations  of  motion  and  the  boundary  con- 
ditions, in  a unified  approach  for  a single  layer  plate. 

A formulation  similar  to  Habip' s for  a general  nonlinear 
sandwich  plate  problem  is  non-existent  in  the  literature. 
Therefore,  it  seems  both  timely  and  important  to  make  an 
attempt  to  fill  this  void.  The  analysis  presented  in  the 
first  part  of  this  dissertation  (Chapter  2 ) has  been  de- 
voted to  making  such  an  attempt. 

Chapter  3 contains  a derivation  of  the  buckling 
equations  by  simplifying  the  results  obtained  in  Chapters 
2 and  3.  These  simplified  equations  have  been  obtained 
by  Kim  [42]  who  has  considered  the  problems  of  combined 
axial  and  transverse  loads,  and  buckling  of  a simply  sup- 
ported rectangular  sandwich  plate.  Thus,  it  has  been  shown 
that  Kim's  equations  can  be  obtained  as  a special  case  of 
the  general  problem  treated  in  Chapters  2 and  3. 

Though  the  elastic  bending  and  elastic  stability 
problems  of  sandwich  plates  have  received  the  attention  of 


numerous  investigators,  the  viscoelastic  properties  of 
the  materials  have  been  neglected.  As  mentioned  in  the 
first  section  of  this  chapter,  the  facings  are  often  made 
of  reinforced  plastic  and  the  core  material  can  be  an 
expanded  material  such  as  cellulose  acetate,  synthetic 
rubber,  etc.  All  these  materials  exhibit  viscoelastic 
behavior  which  may  considerably  influence  the  bending  and 
stability  of  sandwich  plates.  The  author  could  find  very 
few  papers  on  the  viscoelastic  behavior  of  sandwich  plates, 
and  almost  all  of  them  are  very  recent.  This  suggests 
that  the  importance  of  this  aspect  of  sandwich  plates  has 

just  begun  to  be  properly  realized. 

Chang  and  Fang  [43]  made  a valuable  con t r i bu c i on 

by  investigating  the  flexural  motion  of  a rectangular 
sandwich  plate  with  elastic  facings  and  a viscoelastic 
core.  The  dynamic  response  to  a unit  step  function  was 
calculated  first  and  then  the  general  transient  and  steady- 
state  responses  to  an  arbitrary  input  function  were  ex- 
pressed in  terms  of  the  unit  response  by  means  of  convo- 
lution integral.  Malmstrom  [44]  and  Tang  [45]  analyzed 
the  problems  of  creep  bending  of  sandwich  plates  and  two 
layer  plates  respectively,  but  their  theories  are  of  plas- 
ticity. Ashton  and  Moavenzadeh  [46]  considered  a three- 
layer  viscoelastic  half  space  and  analyzed  the  stress  dis- 
tribution in  each  layer.  Ho w ever,  the  analysis  given  by 


Ashton  is  too  general  (perhaps  because  of  his  interest 
in  Civil  Engineering)  to  be  of  any  significance  in  sand- 
wich plate  theory.  Also,  none  of  the  authors  mentioned 
in  this  paragraph  considered  the  problem  of  instability. 

The  only  paper  that  gives  an  analysis  of  bending 
and  buckling  is  by  Lisy  [47]  who  considered  the  problem 
of  a rectangular  simply  supported  viscoelastic  sandwich 
plate  subjected  to  a uniaxial  compressive  load  and  a 
uniformly  distributed  transverse  load.  However,  several 
assumptions  have  been  made  by  Lisy  [47].  These  are: 

( a ) Only  the  core  is  viscoelastic;  the  faces 
are  elastic. 

(b)  Transverse  strain  ir,  the  core  and,  hence, 
the  possibilities  of  the  symmetric  deflection 
and  buckling  are  neglected. 

(c)  Shear  strains  in  the  facings  are  very  small 
compared  to  those  in  the  core. 

(d)  The  material  of  the  core  exhibits  visco- 
elastic behavior  governed  by  linear  Maxwell- 
Thorn  s o n differential  equation. 

(e)  The  deflection  of  the  plate  is  a single  sine 
function  with  a half  wave  length  equal  to 
the  length  of  the  plate  in  each  x and  y 

d i r e c t i o n . 


In  addition  to  the  above  limitations,  the  critical 
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value  of  the  compressive  load  for  buckling  was  calculated 
by  Lisy  only  at  time  = infinity.  Thus  the  applicability 
of  the  result  obtained  by  Lisy  is  severely  limited  in 
scope . 

In  this  dissertation,  an  attempt  has  been  made  to 
analyze  a similarly  loaded  sandwich  plate  but  without 
making  any  of  the  assumptions  listed  above,  except  assump- 
tion (d).  Linearity  in  the  material  response  function  has 
been  retained  as  the  first  order  approximation  to  the  actual 
case,  while  a better  and  more  complicated  model  will  per- 
haps be  the  subject  of  future  investigations. 

Chapter  4 gives  the  analysis  of  viscoelastic  sand- 
wich plates  derived  from  elastic  sandwich  plates  using  the 
correspondence  principle,  and  Chapter  5 gives  the  numerical 
results  and  conclusions. 


CHAPTER  2 


GENERAL  NONLINEAR  ANALYSIS  OF  ELASTIC  PLATES 


Notation  and  Coordinates 

Throughout  this  dissertation,  Cartesian  coordinates 
have  been  used.  In  most  places,  especially  in  the  first 
part  of  this  dissertation  where  long  expressions  are  in- 
volved, indicia!  notation  has  been  used  for  the  purposes 
of  compactness  and  convenience.  If  the  components  of  a 
vector  a in  the  x,  y,  z coordinates  are  ax,  a , a^  re- 
spectively, in  indie i a 1 notation  these  will  be  denoted  by 
a ^ , i = 1,  ^ , 3.  Thus,  a ^ — a -|  , ^y  — ^ 2 3 ^ ^ ~ ^3' 

Also,  the  usual  summation  convention  lias  been  used  in  con- 
nection with  indicial  notation.  According  to  this  conven- 
tion, whenever  the  same  letter  subscript  occurs  twice  in  a 
term,  that  subscript  is  to  be  given  all  possible  values  and 
the  results  are  to  be  added  together.  The  symbol  E may 
then  be  omitted  because  the  repeated  subscript  indicates 
the  summation.  For  example,  in  three  dimensions, 

3 i 2 b i = a 1 2 b 1 + a 22  b2  + a32b3 
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The  Latin  indices  (i,  j,  k,  1,  etc.)  assume  the 

values  1,  2 , and  3,  while  the  Greek  indices  (a,  3 , Y , 6 , 

etc.)  assume  the  values  1 and  2 only.  Thus,  these  two  sets 

of  indices  can  be  used  for  three-dimensional  quantities 

and  two-dimensional  quantities  respectively.  In  addition, 

the  usual  notation  A . . which  denotes  9 A - has  been  used 

^5  J ]_ 

d X 

throughout  this  dissertation.  j 

In  the  present  analysis,  the  Modified  Hellinger- 
Reissner  Principle  has  been  taken  as  the  starting  point 
from  which  subsequent  analyses  leading  to  various  special 
cases  have  emerged  in  a number  of  steps.  The  Hellinger- 
Reissner  Principle  and  the  Principle  of  Minimum  Potential 
Energy  can  be  interpreted  as  special  cases  of  the  General- 
ized Principle  of  Minimum  Potential  Energy  [48].  The 
He! 1 i nger-Rei ssner  Principle  as  given  in  [48]  involves 
linear  strain-displacement  relations,  while  the  Modified 
Hell i nger-Rei ssner  Principle,  used  by  Habip  [41]  on  the 
basis  of  [49]  is  also  applicable  to  non-linear  strain  dis- 
placement relations. 

Three-Dimensional  An  a 1 y s i s 

Modified  H e 1 1 i n q e r-Reissner  Prin c i p 1 e 
If  for  an  elastic  body 

0V,  dV  = Volume  of  the  undeformed  body,  related 
element  of  volume 


1G 


° A s ( t ) ’ 


dA 


oAV(t) ’ dA 


V . 


po 
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fi 
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F . 


£* 


s 
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= Total  surface  area  of  the  undeformed 
body  over  which  stresses  are  pre- 
scribed, related  element  of  area 
= Total  surface  area  of  the  undeformed 
body  over  which  displacements  are 
prescribed,  related  element  of  area 
= Components  of  the  displacement  vector 
referred  to  the  base  vectors  of  the 
undeformed  body 
= Strain  tensor 

= 2nd  Piola-Kirchhoff  stress  tensor  as 
defined  in  [50] 

= Density  of  the  undeformed  body 
= Components  of  acceleration  vector 
referred  to  the  base  vectors  of  the 
undeformed  body 

= Components  of  the  body  forces  per 
unit  mass  of  the  undeformed  body  re- 
ferred to  the  base  vectors  of  the 
undeformed  body 

= Strain  energy  function  per  unit  volume 
of  the  undeformed  body 
= Components  of  the  stress  vector  0S 
referred  to  the  base  vectors  in  the 


undeformed  body 


(1) 


I / 


6 V = Variation  of  a quantity  V 

And  if  tildes  (%)  indicate  prescribed  quantities,  then 
the  Modified  Hellinger-Reissner  Principle  as  given  in  [41] 
asserts  that  the  variational  principle 


(•V-F^VidV  - dffejYjj-e*)  *-(iSij(»1tjt»j,1+»k.iVk.j)dV 


°v 


+ 1 Si*VidA  + 


si*Cvi-Vi)<la 


(2) 


oA 


s(t) 


4V(t) 


where  y..,  S..,  V.  and  S.*  are  varied  independently,  is 
i j i J i i 

equivalent  to  Cauchy's  first  law  in  0V,  to  the  stress  bound- 
ary condition  on  the  part  0As^^  of  the  boundary,  to  the 
displacement  boundary  condition  on  the  remaining  part  0Ay, 
and  to  the  stress-strain  and  strain-displacement  relations 
in  o V,  when  the  symmetries  of  y..  and  S..  are  both  used. 

I J 1 J 

Habip  [41]  has  demonstrated  the  proof  of  the  above 
stated  theorem.  In  this  dissertation,  the  following  strain 
displacement  relations 


a{3 


y 


a 3 


- + Va  + V3,aV3,el 

* *<Va,3  + V3,al 
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(3) 


have  been  used  instead  of  the  exact  relations 


. + V . . 

> J J , i 


(4) 


Relations  (3)  have  been  used  by  many  authors,  e.g.,  [19], 
[51],  and  are  realistic  considering  the  degree  of  non- 
linearity and  therefore  complexity  involved  for  the  corres- 
ponding gain  in  accuracy.  However,  for  very  large  defor- 
mations, Equation  (4)  must  be  used. 

Corresponding  to  the  strain-displacement  relations 
adopted  here  and  given  by  (3),  the  Modified  Hellinger- 
Reissner  Principle  given  by  (2)  changes  to 
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s Jv  » + v, 

a3V  a, 3 3, a 


) dv  - 
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S33V3,3dV 
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oV 


J 

O V 


oAs(t) 


(5) 


Av(t) 


Equation  (5)  is  the  governing  equation  of  the 
present  analysis.  In  both  sides  of  Equation  (5),  the 
operations  of  variation  (6)  are  carried  out  within  the 
integrals  and  Green's  theorem  is  made  use  of  to  convert 
volume  integrals  into  surface  integrals.  To  avoid  writing 
down  very  lengthy  expressions  unnecessarily,  only  a sam- 
ple calculation  for  one  of  the  expressions  in  Equation  (4) 
is  shown  below,  calculations  for  the  other  expressions 
being  quite  similar. 
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(Ga36Va),3dv  - [fSaBl.B6VadV 
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aB 


;nD6V  dA  + ( S 0OnQ6V  dA  + S .0n,6V  dA  + S ^cn,6V  dA 
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oAs(t)  oAv(t)  oAs(t)  oAv(t) 


(contd. ) 
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a3  a 

Tv 

where 
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-/< 


(S  -6V  ) dV  - (S  _)  5 V dV 


,V 


aB  a 


(6) 


oA 


(t) 


° n i 


= Total  surface  of  the  body 

= Cosine  of  the  angle  between  the  out- 
ward normal  to  the  surface  of  the 
body  at  the  point  concerned  and  the 
coordinate  axis  referred  to  the 
undeforned  body . ( 7 ) 

After  carrying  out  the  calculations  for  all  the 
expressions  in  Equation  (5)  in  a way  similar  to  that  shown 
in  (6)  and  collecting  all  the  expressions  on  one  side, 
Equation  (5)  takes  the  following  form: 


(A)  + (B)  + (C)  = 0 

where 


(8) 
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(B) 
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(9) 


Up  to  this  point,  this  analysis  has  been  very  gen- 
eral and  three-dimensional  with  no  reference  to  a plate- 
like two-dimensional  problem.  Equation  (8)  with  the  aid 
of  (9)  can  be  applied  to  any  three-dimensional  elasticity 
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problem.  To  do  so,  one  would  have  to  use  the  fundamental 
lemma  of  the  Calculus  of  Variations  [52]  and  set  each  in- 
dividual integrand  contained  in  Equation  (8)  equal  to  zero. 
A complete  set  of  equilibrium  equations,  elastic  equations, 
compatibility  relations  and  boundary  conditions  would  be 
obtained.  For  two-dimensional  problems,  this  has  been 
done  in  the  subsequent  chapters,  and  the  procedure  for 
the  three-dimensional  problem  is  similar  to  the  two-dimen- 
sional one.  Since  our  aim  is  to  analyze  plate  problems, 
from  this  point  on,  our  analysis  would  be  two-dimensional. 

Tw o - Pi  mens  i ona.l  Analysis 

Figure  (1)  shows  a body  having  two  parallel  and 
flat  surfaces,  and  surrounded  by  a surface  perpendicular 
to  the  two  flat  faces.  The  thickness  is  small  compared 
to  other  dimensions  of  the  body,  so  we  can  consider  it 
to  be  a plate.  The  plate  can  have  many  layers  or  just  a 
single  layer. 

The  coordinate  axes  x ^ , x^ , and  x^  are  taken  such 
that  the  x ^ x ^ plane  coincides  with  the  middle  plane  of  the 
plate. 

Various  quantities  are  defined  as  follows: 

(A  subscript  "0"  refers  to  the  undeformed  state.) 

0V  = Total  volume  of  the  plate 

0A  = Total  area  of  the  middle  surface 
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Figure  1 . 


Plate  with  arbitrarily  specified  loads  a r, d 
displacements  over  its  surface 
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where 


Total  surface  of  the  plate 
oAs(t)  + oAV(t) 


°As(t)  anc*  oAv(t)  *iave  been  defined  in  (1) 
oAs(t)  = ° As ( U+L ) + oAcs 

°As(U+L)  = Total  portions  of  the  upper  and  lower 

surfaces  of  the  plate  over  which  stresses 
are  prescri bed . 

= oAs(U)  + ° As ( L ) 


o 


A 


s ( U ) 


° As ( L ) 


o 


A 


cs 


oAV(t) 

oAV(U+L) 


Portion  of  the  upper  surface  over  which 
stresses  are  prescribed 
Portion  of  the  lower  surface  over  which 
stresses  are  prescribed 

Portion  of  the  edge  surface  of  the  plate 
over  which  stresses  are  prescribed 

oAV(U+L)  + oAcV 

Total  portions  of  the  upper  and  lower 
surfaces  over  which  displacements  are 
prescribed 

°AV(u)  + oAv(L) 


oAV(U) 

oAV(L) 


= Portion  of  the  upper  surface  over  which 
displacements  are  prescribed 
= Portion  of  the  Lower  surface  over  which 
displacements  are  prescribed 
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0Acy  = Portion  of  the  edge-surface  over  which 
displacements  are  prescribed. 


In  our  case, 


°As ( U)  oAs(L) 

in  magnitude 

oAs 

(say) 

Also, 

oAV(U)  ° Av ( L ) 

in  magnitude 

° A v 

(say) 

(10) 


Now  the  results  of  Equation  (8)  can  be  applied  to 
the  plate  considered  above.  Due  to  difficulties  in  carry- 
ing out  very  lengthy  calculations  for  different  expressions 
all  t h rough  simultaneously,  the  calculations  resulting  from 
the  expressions  (A),  (B),  and  (C)  of  Equation  (8)  would  be 
carried  out  separately.  Expression  (B)  would  be  considered 
first. 


In  the  expression  for  (B)  given  in  Equation  (9), 
each  of  the  three  volume  integrals  is  converted  into  a 
double  integral  consisting  of  an  area  integral  over  the 
middle  surface  of  the  plate  and  a simple  integral  through 
the  whole  thickness  of  the  plate.  The  surface  integrals 
like  0As£t)  are  broken  up  into  two  parts,  °As(y+[_)  anc* 
°^cs'  The  sur^ace  integrals  associated  with  0Acs  are  then 
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converted  into  a double  integral,  one  of  which  is  a line 
integral  along  0 C s , the  edge  of  the  middle  plane  over 
which  stress  resultants  are  given,  and  the  other  one  is 
a simple  integral  through  the  whole  thickness  of  the  plate. 
Then,  noting  that 


and 


on^  = 0 on  the  upper  and  lower  surfaces 
on^  = 0 on  the  edge, 


(B)  can  be  expressed  as  follows: 


(B)  = (D)  + (E)  + (F)  + (G) 
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(ds  = Element  of  length  along  DCs  or  DCy) 


(12) 
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The  results  obtained  so  far  are  applicable  to 
plates  consisting  of  any  number  of  layers  having  differ- 
ent properties  and  geometry.  But  from  here  on,  the  analy- 
sis would  depend  on  the  type  of  plate  considered  viz.  sin- 
gle layer  or  multiplayer  plates,  thin  or  moderately  thick 
plates,  etc.,  and  also  on  the  type  of  displacement  func- 
tions Va  and  Vg.  Habip  [41]  has  demonstrated  the  applica- 
bility of  these  results  to  a single  layer  plate. 

In  this  dissertation,  a three-layered  plate,  popu- 
larly known  as  a sandwich  plate  will  be  considered. 

Sandwich  Plate 

A sandwich  plate  consisting  of  three  layers  is 
shown  in  Figure  (2).  The  face  layers  are  much  thinner 
than  the  core.  All  layers  are  uniformly  thick  throughout. 
The  two  faces  are  of  the  same  thickness  t ^ , the  core  thick- 
ness being  t . 

a c 

The  present  analysis  is  capable  of  treating  mixed 
boundary  value  problems.  On  the  upper  surface,  the  line 
AB  separates  the  two  regions  0A^u^  and  0As^  over  which 
displacements  and  stresses  respectively  are  prescribed. 

On  the  lower  surface,  the  line  CD  separates  the  correspond- 
ing regions.  The  line  CD  is  located  exactly  below  the  line 
AB.  Also,  at  any  two  points  located  on  the  upper  and  lower 
surfaces  of  the  plate  and  having  the  same  x-j  and  x ^ 
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Figure  2.  Sandwich  plate  of  a general  shape 
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coordinates,  the  same  type  of  quantities  (stresses  or  dis- 
placements) are  specified.  Therefore  0Ay^  is  equal  to 
°%(L)  both  in  size  and  shape,  and  the  former  is  located 
exactly  above  the  latter,  as  shown  in  the  figure.  The 
same  thing  holds  for  0A^yj  and  0As^.  Finally,  the 
lines  EF  and  GH  which  separate  the  portions  0 A c s and  0 A c y 
of  the  edge  surface  of  the  plate,  are  straight  lines  per- 
pendicular to  the  middle  surface  of  the  core.  These  con- 
ditions are  necessary  because  they  assure  that  the  type 
of  prescribed  quantities  (stresses  or  displacements)  do  not 
change  across  the  thickness. 

In  addition  to  the  conditions  stated  above,  the 
following  assumptions  are  made: 

Al)  The  transverse  normal  strain  in  the  face  layers 
is  negligible  compared  to  that  of  the  core. 

A2 ) Displacements  in  each  layer  are  linear  func- 
tions of  distances  from  the  median  plane  of 
that  layer. 

A3)  The  core  is  attached  to  the  faces  securely. 

For  almost  all  the  common  shapes  and  material  prop- 
erties of  a sandwich  plate,  assumption  (Al)  is  valid. 

Typical  geometry  and  material  properties  of  commercially 
used  sandwich  plates  can  be  found  in  [1]  and  [53].  Assump- 
tion ( A2 ) has  been  used  by  [10],  [11],  [12]  and  many  others. 
This  is  an  approximation  to  a power  series  containing 
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an  infinite  number  of  terms  for  the  displacement  functions. 
The  validity  of  this  assumption  depends  on  the  amount  of 
deformation  produced  and  the  accuracy  needed.  However, 
retaining  even  a single  nonlinear  term  in  the  expressions 
for  the  displacements  give  rise  to  inordinately  cumbersome 
equations  with  no  appreciable  gain  in  accuracy.  Thus  this 
assumption  is  justified.  Assumption  (A3)  rises  out  of  the 
necessity  of  continuity  of  displacements  at  the  interfaces 
of  different  layers. 

Three  sets  of  coordinates  as  shown  in  Figure  (2) 
are  chosen  such  that  their  origins  lie  on  the  same  straight 
line  perpendicular  to  the  u n deformed  middle  plane  of  the 
core.  The  planes  containing  (x^‘,  x2  1 ) , (x-j,  x^)  and 
(x] " , x2" ) axes  coincide  with  the  middle  planes  of  the 
upper  layer,  the  core  and  the  lower  layer  respectively, 
and  the  x ^ 1 , x^  and  x^"  axes  are  on  the  same  straight  line 
and  in  the  same  direction. 

The  single-primed  and  double-primed  quantities 
refer  to  the  upper  and  lower  faces  respectively.  The  un- 
primed quantities  refer  to  the  core.  Then,  the  following 
relations  hold: 

x3 ' = x3  - A(tc  + tf) 
x3"  = x3  + i(tc  + tf) 


X 


a 


(13) 
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Assumptions  (Al)  and  (A2)  lead  to  the  following 
relations  between  different  displacements: 


V = U'a  <V  + X3  ' W 
V3'  = V'V 

V ’ V<XB>  + VV  <XB> 


V 


U3” *XB ' 


a 


W + X3^a  (XB> 


V3  = U3(xe)  + X3H'3(X6) 


(14) 


where  U U0' , T ' , U 11  , U ",  T"  , U , U,.  T and  'f , are  called 

a ' 3 a a o a a 3 ' a 3 

displacement  functions  and  are  functions  of  x-j  and  x ^ only. 
There  are  16  displacement  functions  in  Equation  (14).  But 
the  assumption  (A3)  gives  rise  to  the  following  relations 
among  these  quantities: 
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X 1 3 X 2 5 

X 
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if)  = 
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Vx  !■ 

x 2 5 

X3 

t c \ 
2 1 

( X 1 9 X 2 5 

X3  ' = 

_li)  = 
2 

v3(xl , 

x2  , 

x3 

tc  X 
2 

( X | 9 X 2 5 

x3"  = 

— ) = 
2’ 

Vv 

x 2 5 

X3 

tc  X 
2 

( X -|  5 X 2 9 

x3"  = 

— ) - 
2 ' 

v3(xl , 

x 2 ? 

x3 

tc  X 
2' 

(15) 

There  are  six  equations  in  (15).  Thus  the  number  of  in- 
dependent displacement  functions  is  reduced  from  16  to  10 
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by  combining  (14)  and  (15).  Application  of  conditions 
(15)  to  (14)  gives 


a 


= |~(U'  -U  ) 

t ^ a a 2 a 

tc 

U3  + '2  ^ 3 


II 


¥ 


a 


U 


3 


(U  -U 

a a 


(16) 


Substitution  of  (13)  and  (16)  into  (14)  gives  the  expres- 
sions for  the  displacements  in  terms  of  the  displacement 
functions  as  follows: 


2 r_.  1 /j_  . J.  \1  n T,  2 I-..  1 /j.  , J.  \1  C 


a 


Uk  t1  + 17  (x3-*(tc+tf»  1 - UcFf  tx3-4(tc+tf»-  Vft  fx3--(tc+tf» 

I j + . _~\|/ 

u3  2 3 


V " 

a 


C’-nfs  +*(tc+tf)} ] + u0t~{ x3+f{tc+tf)}  -Vr:{Vi<tc+tf>} 


a 


V3" 


u3  - -ih 


V = U + x,T 
a a da 


V3  U3  + x3T3 


(17) 


The  expressions  given  by  (17)  are  then  substituted 
into  the  expressions  (D),  (E),  (F),  and  (G)  given  by  (12) 
and  in  the  expression  (C)  given  by  (9).  The  integrations 
are  then  carried  out  through  the  thickness  of  the  plate. 
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The  expressions  that  emerge  from  each  of  the  quanti- 
ties (D),  (E),  (F),  etc.  are  excessively  long  in  the  inter- 
mediate steps  making  it  extremely  inconvenient  to  carry  out 
calculations  of  all  of  them  simultaneously.  Therefore, 
calculations  would  be  shown  for  (D)  only,  those  for  the 
rest  being  similar. 


like  (D)  alone,  involve  very  long  expressions  in  the  inter- 
mediate steps,  and  require  a great  amount  of  writing  space 
to  show  the  procedure  followed.  On  the  other  hand,  omit- 
ting all  the  intermediate  steps  may  make  it  very  hard  to 
know  hew  the  final  results  were  determined.  Therefore 
only  those  intermediate  steps  which  are  essential  to  keep 
track  of  the  procedure  followed  will  be  shown. 


and  using  definitions  (10),  the  first  integral  of  the  ex- 
pression (D)  is  given  as: 


But  the  calculations  for  each  of  these  quantities, 


Noting  that 


on  the  upper  surface 


1 on  the  lower  surface 
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{ 2 ( S ' * ) , 


2(S'3)h}6U;dA+r{2(s^)_h  + 2(s^3)_t,}6U;dA 


% 


a*  h 


,A_ 


'Xj 


'Xj 


{(S^>-h+'Sa3>-h+<S;*)h-(Sa3>h>Syadfl 


oA. 


+J<T  ^*>-h  - \ <S'a*>h  + T <S:3)-h  + T-  ^S^V*  (,8) 

oAs 

where  the  following  notation  has  been  used: 

(x)  = value  of  x at  x,  = y. 

y ^ 

In  a way  similar  to  the  one  just  shown,  the  second  integral 

in  the  expression  of  D is 


{S3*_'S36V3,3+S33^°n3}oV3dA  f{^S3*^-h  + 'S3*^  + ^S30^-h  U3,f 


>A 


s(U+L) 


2^S3B^-h  T3,3+  ^S33^-h  “ ^S33^hU3,f 


2^S3B^h  ¥3,b 


( S33 ) h ^ 6U3dA 


+ f {-1(S3*)-h  + 2 ^S3*^h  2 ^ S3 g ^ - h U3,B+  4 ^V-h  ^3,3 


1 


2 ^ S33 ^ -h  " 2 ^S33^U3,3  4 ' ^33^3, 3"  2 ^S33^h}  SV3dA  (19) 
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The  third  integral  in  the  expression  for  (D)  is 
considered  next.  The  integral  through  the  whole  thick- 
ness of  the  composite  plate  between  the  limits  -h  and 

+h  is  divided  into  three  integrals,  the  first  one  b e - 

tc 

tween  the  limits  -h  and  - — ^ , the  second  one  between 
tc  ^ 

the  limits  - — ~ and  — and  the  last  one  between  the  limits 

tc 

— 2 and  h,  corresponding  to  the  lower  face,  the  core  and 
the  upper  face  respectively.  For  each  layer,  the  appro- 
priate expressions  for  displacements  are  substituted  in 
the  integral  for  that  particular  layer.  Thus, 
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(20) 
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One  important  point  regarding  notation  is  to  be 
noted  in  the  above  equation.  The  left  hand  side  is  a 
general  expression  which  i s a p p 1 i cabl e to  a single  layer 
or  a multilayer  plate.  So  all  quantities  in  the  left 
hand  side  are  unprimed.  But  the  right  hand  side  is  for  a 
three-layered  sandwich  plate,  so  the  unprimed  quantities 
in  the  right  hand  side  refer  to  the  core,  as  stated  earlier 
Now,  the  following  changes  of  variables  and  limits 
of  integration  can  be  made: 


jc 
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1 
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Y(xo)  dx3 


y<xJ 


dx3" 


-h 


1 
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Y<xa>  dx3 
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y 


2 

Y ( x ) dx. 


a 


(21) 


If 

2 


[using  Equation  (13)] 


Then,  substituting  Equation  (17)  into  (20)  and 
using  relations  (21),  the  following  results  are  obtained: 


[The  3rd  integral  in  the  expression  of  (D)] 

l(S;B.6+P”"("Fa-“fa)H6Ua(,-f7V)  + V^alT  V1  dX’"dA 
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(22) 
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Similar  procedure  is  applied  to  the  fourth  integral  in 
the  expression  of  (D).  The  last  integral  in  the  expres- 
sion of  (D)  is  now  considered: 
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(23) 


Then  the  final  expressions  given  by  (18),  (19),  (22),  (23) 
and  the  one  for  the  fourth  integral  in  the  expression  of 
(D)  are  substituted  into  the  first  equation  of  (12).  After 
that  the  remaining  integrations  are  carried  out  and  the 
terms  within  the  same  integral  and  associated  with  variation 
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of  the  same  quantity  are  grouped  together. 

1ike  ‘ T7  <Si3>^c  and  17  (Sa3>tc  are  made 

2 f 2 

other  because  at  the  interface  between  the 

the  core  (S^)^.  = <So3)t  . 

2 2 


Then  the  terms 
to  cancel  each 

upper  face  and 


The  next  step  in  the  present  analysis  is  to  define 
a number  of  quantities  as  follows: 
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Using  definitions  (24),  the  final  expression  for  the 
quantity  (D)  becomes: 
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(25) 


The  quantity  (.F)  given  in  (12)  is  calculated  in  a 
way  exactly  similar  to  that  for  the  quantity  (D),  and  the 
resulting  expression  for  (F)  is  exactly  the  same  as  for 
(D)  given  in  (25)  except  that  the  quantities  and  S^* 

are  to  be  replaced  by  and  S^*  respectively. 

The  procedure  of  evaluating  the  quantity  (E)  is 
similar  to  that  for  (D)  just  shown.  Only  the  final  result 
is  given  below: 
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where,  the  following  definitions  have  been  used  in  addition 
to  those  given  in  (24): 
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The  expression  (G)  given  in  (12)  is  evaluated  in 
a way  which  is  similar  to  that  for  (E)  just  shown  above. 
Also,  the  resulting  expression  for  (G)  is  exactly  the  same 
as  that  for  (E)  given  in  (26)  except  that  the  tildes  (^ ) 
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are  removed  from  all  quantities  in  which  they  occur,  in- 
cluding the  definitions  (27). 

The  quantity  (C)  given  in  (9)  can  be  evaluated  by 
the  procedure  used  already  for  the  quantities  (D),  (E), 
etc.,  and  the  resulting  expression  for  (C)  is  given  below: 
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Finally,  the  quantity  (A)  given  in  (9)  is  now  con- 
sidered. Like  the  procedure  followed  previously,  each  sin- 
gle volume  integral  is  first  converted  into  a double  inte- 
gral, one  of  which  is  an  area  integral  and  the  other  one  of 
which  is  an  integral  through  the  thickness.  The  procedure 
is  demonstrated  with  the  first  integral  in  the  expression 
of  (A)  given  in  (9). 
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The  following  forms  are  assumed  for  the  strains 
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n = 0 


where  M,  M1  and  M"  are  positive  integers.  The  quantities 

v „y 1 n and  „y''  are  obtained  in  terms  of  the  displace- 

n "aB  n’aB  n'aB 

ment  functions  given  in  (14)  after  (30)  is  substituted  into 
(29)  and  into  the  second  and  the  third  integrals  in  the  ex- 
pression of  (A). 
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The  fourth,  the  fifth  and  the  sixth  integrals  in 
the  expression  of  (A)  give  the  relationships  between  vari- 
ous stress  resultants  (N  „ M _etc.)  and  the  strain  func- 
tions  (nYag)  through  the  strain  energy  function  per  unit 
area  (e)  for  each  layer.  The  procedure  is  demonstrated 
with  the  fourth  integral. 


9c*  \ . 

9V  ap 


9e*'\  x ■■ 

Kb  Y“6  3 


h 

+ ^Sc*B  §77^  °YaBdX3  + J ^Sa3  ” WV  6YaB  dX3^  dA 

a 3 ^ j.3 


tc 

2 


:3i) 


7 

2 


"«-*  I Kb  - ig>  s''ae  dx3  -J  Kb  - !*£’  dV 

-h  - tf 


(32) 


[by  applying  (21 )] 

From  (30),  the  following  result  is  obtained: 


9c*" 


9c*" 

9 Y " 

1 y v 


9Y"yy 

><b 


9c*" 

9y " , 
1 y v 


Kv 


+ n Y " ) 

1 yv 


9e*  " 


(33) 


Similarly, 


H / 


3e*"  _ „ 3e*" 

VJd  ' X3  Ke 


(34) 


Also, 


5 <6 


VaB 


(35) 


With  the  aid  of  (33),  (34)  and  (35),  the  right  hand  side 
of  Equation  (32)  takes  the  form: 


_ *n 

-Yrr~) 

° o.& 


(x,"  S" 


3e*"  ^ 


Va 


. X 

" ' °l'a.B 


] dx, 


o r 


-h 


(S)' 


a 8 


Pic*" 

~n— ) 6y"0  dXo 

3YaS  “S  3 


(N‘‘ 

a8 


6<B  + (M«B 


3e" 

3,  y r 

i a8 


] VaS 

(36) 


(37) 
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The  second  and  third  integrals  on  the  right  hand 
side  of  Equation  (31)  are  evaluated  in  a similar  way.  The 
procedure  is  then  repeated  for  the  fifth  and  the  sixth  in- 
tegrals in  the  expression  of  (A)  given  in  (9). 

Now  the  resulting  expressions  of  the  quantities 
(D),  (E),  (F)  and  (G)  are  added  to  obtain  the  quantity 
(B),  and  then  to  the  quantity  (B)  the  values  of  the  quan- 
tities (A)  and  (C)  are  added.  This  combined  expression  is 
equated  to  zero  as  per  Equation  (8).  Since  all  the  inde- 
pendent quantities  can  be  varied  independently  of  each 
other,  then  by  the  fundamental  lemma  of  Calculus  of  Varia- 
tions, each  expression  within  a separate  integral  and  asso- 
ciated with  the  variation  of  an  independent  quantity  must 
be  equal  to  zero.  All  the  governing  equations  for  a sand- 
wich plate  are  then  generated. 

At  this  point  it  is  to  be  noted  that  the  original 
three-dimensional  sandwich  plate  problem  has  now  been  re- 
duced to  a purely  two-dimensional  problem  [Figure  (3)]. 

Figure  (3)  shows  a flat  surface  of  total  area  0A 
and  bounded  by  a closed  curve  QC.  The  total  area  CA  is 
divided  into  two  parts,  0Ay  and  0AS , by  a line  AB  (which 
can  be  a closed  curve,  or  a combination  of  closed  curves  in 
the  general  case),  and  over  these  two  parts,  displacements 
and  stresses,  respectively,  are  prescribed.  The  periphery 
0C  is  divided  into  two  parts,  DAC  and  CBD,  which  are 
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Figure  3 . 


Two-dimensional 


mixed  boundary  value  problem 


called  QCs  and  0 C y respectively.  Prescribed  quantities 
are  stresses  over  °Cs  and  displacements  over  0CV,.  It  is 
therefore  a mixed  boundary  value  problem. 


Equilibrium  Equations  Over  0 A 


a£,B 


- v * F;-f; + 4 <«> + 24a 


= o 


N"Q  d - a - Q")  + FM-f"  - ) + 2(S'J  . = 0 

a3»B  t^  aB,B  a a a t^  a a a*  -n 


4 <«S,S  - Qa  + W + "aB.B  + Fa  ' fa 

- 4 <«-«>  - ■ ° 
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~ (Q'-M1  + Q"-M"  ) + M - 0 + M 

t rr  ' XX  00,8  XX  00,0'  00,0  XX  ( 


m 

a a 


- (M,-m'+M"-m")  (S"J  . } -0 

t^  a a a a 2 a*  h ' a* -h 


tc 


{N'  (U,  + -4  J + N"  (U 


a0v  3,0  2 3,0 


d)  + N Q U , o + M . 


a0  v 3,0  2 '3,0'  aB  3,0  00*3,0 


a. 


+Q 1 +Q"+Q  } + FJ,-f ‘ +F~-f ';+Fo-fo+  (S'  ).  + (S"  ) . = 0 

a a Ha  ,a  333333'  3*  h 3*  -h 


and 


%No0^U3,0+  2 ¥3,0^  + Qo},o  " 2 No0(U3,e"  2 'F3,0)  + Qa}  ,a 


tc 


+ (M  0IL  „+K  >f'  +T  ) - N„  + (F^-f  ‘ -F"+f") 

a0  3,0  a0  3,0  a ,a  33 


3 3 


+ M3  m3  + 2 (S3*)h  (s3*)_h  0 


(38) 


Equilibrium  equations  over  0 A y are  exactly  the  same 
as  those  on  0As  except  that  the  tildes  (a.)  are  removed  from 
all  the  quantities  in  which  they  appear. 


Boundary  Conditions 


(1)  On  0C 


(S’),  - ~ (t1),  - (N ' - ) 0n„  - 0 

ah  t^  ah  a0  t^  a0  0 


^ p % p 

(S")  . + t~  (t")  , - (N"q  + t~  M"0 ) on0  - 0 
' a -h  tf  a -h  a0  tf  a0  8 


4 «\l>h  - h-u}  + X-  Tf  <Ma6  • V »"e  - V"e  ■ 0 


^{(*d>h  + '‘d>-h  } -X-rf  + •&>  °nB  + V"b  - 0 


-V 


% 'b 


, tr  tr 

^ Ml  '-Mil 


<Si>h  + S3  +<S3>-h  - ‘^bW  u3,8  ' <“?  V T NdB 


+ M J H'o  o + Q1  + Q + Q"  } on  = 0 
a$  3,3  tx  a a a 


and 


tc  ^ 


{ ( S3 ) h + ^SP-h  } + t3  “ ^ 2 { Na8^U3,3  " f '3,3 ^ + Q«v  +Mo;BU3,3 


+V3,3  + Ta  + T{Na6(U3,e+-T’U3,3)+Qa}]ona=  ° (39) 


(2)  On  0C v: 


Boundary  conditions  on  o0y  are  exactly  the  same  as 
those  on  QCs  except  that  the  tildes  (^)  are  removed  from 
all  the  quantities  in  which  they  appear.  In  addition,  the 
following  conditions  hold  on  oC^: 


u1  = 

a 

% 

U'  , 

a 

U = 
a 

% 

U , 

a 

U" 

a 

= U", 

a 

U3 

r\j 

U3’ 

¥ 

a 

<? 

II 

(3) 

On 

0 Ay  : 

U'  - 

a 

'b 

U1  , 

a 

U = 
a 

'b 

U , 

a 

U" 

a 

= U", 

a 

U3 

a. 

U3  ’ 

Y 

a 

a 

c5^ 

ii 

(40) 


'b 

'f. 


(41  ) 
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Strain  Functions  in  Terms  of  Displacement  Functions 


°Y 


' = i(u'  Q + u;  + u'  u'  .) 

a3  a , 3 3 , a 3, a 3,3 


rig  = *<^,6  + 


> Y 


n _ i 


a3 


(U"  . + U"  +U"  U"  D) 


a , 


,a  3 ,a  3,| 


v " = A ( w " + w " ) 

1 ' aB  a,3  B,a; 


°ya3  2^Uct,3  + U3,a  + U3,aU3,3^ 


Y o “ £(¥  + Uo  o + U,  ) 

1 aB  a, 3 3 ,a  3, a 3,3  3,3  3 , a 


v = iu/  w 

2 1 a3  2 3 , a x3,3 


oY ’ o - 1 + U ’ ) 

' a3  2 a 3, a 


oY"q  = hiV"  + U"  ) 
a3  a 3 ,a 


oY  q = 4(Y  + U,  ) 
a3  a 3 ,a 


lYa3  21 3, a 


Y . 


,y33  ^ 3 
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and 

n^B  = n^6 


nYc.3  = 0 for  n>2 


nYa3  = 0 for  n > 3 

nYa3  = "Ya3  = nY33  ' 0 forni1 
n Y 3 31  = nY33"  = ® for  values  of  n 


(42) 


Relations  Between  V arious  Stress-Resultants  and  the 
Strai n Functions 
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9c" 
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KaB= 
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3lYaB 

VaB 
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N33 

9e 
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9°Y33 

(43) 


Equations  (33)  through  (43)  constitute  a comp! ete  set 
of  equations  necessary  and  sufficient  to  solve,  in  principle, 
any  mixed  boundary  value  problem  of  sandwich  plates. 

The  equilibrium  Equation  (38)  in  a slightly  more 


r /» 

j 


general  form  (because  of  using  the  exact  Equation  (4) 
instead  of  the  approximate  Equation  (3)  which  has  been 
used  in  the  present  analysis)  were  obtained  by  Ebcioglu 
[20]  using  a different  approach.  The  reason  for  choosing 
the  approximate  relations  (3)  instead  of  the  exact  rela- 
tions (4)  has  been  explained  before.  Also,  the  present 
analysis  could  be  carried  out  using  the  exact  Equation 
(4)  without  any  difficulty. 

The  advantages  of  using  the  present  analysis  over 
the  one  by  Ebcioglu  [20]  are:  (a)  formulation  of  a more 

general  mixed-boundary  value  problem,  and  (b)  a unified 
treatment  of  a general  sandwich  plate  problem.  This  means 
that  the  strain  functions  in  terms  of  the  displacement 
functions  given  by  Equation  (42)  and  the  stress  resultants 
in  terms  of  the  strain  functions  given  by  Equation  (43) 
are  obtained  as  results  of  the  present  analysis  along  with 
the  equilibrium  equations  and  boundary  conditions,  unlike 
the  analysis  given  in  [20].  Thus  the  present  analysis  can 
be  considered  to  be  a more  general  one  than  the  previous 
work  by  Ebc i oglu. 

Constitutive  Equations  and  Stress  Resultants  in  Terms  of 
Displacement  Functions  for  Isotropic  Material 

The  analysis  presented  so  far  does  not  consider  any 
specific  constitutive  equation  (i.e.,  stress-strain 


55 


relationship)  and  is  therefore  valid  for  any  kind  of  elas- 
tic material  (like  non-homogeneous  and  orthotropic).  This 
section  considers  a sandwich  plate  whose  faces  and  cores 
are  made  of  isotropic  and  homogeneous  materials. 

Both  the  materials  are  made  of  the  same  material 
having  shear  modulus  and  Poisson's  ratio  equal  to  G ^ and 

Vf  respectively,  and  the  core  made  of  a different  material 

having  shear  modulus  and  Poisson's  ratio  equal  to  Gc  and 

respectively. 

For  a linearly  isotropic  material,  the  constitutive 

equations  are: 

S • • = Ay.  . <5 . . + 2uv . . (44) 

i 3 'kk  ij  j v ' ' 

where  A and  y are  Lame-constants,  and 
2vG 

A - 1 - 2 v 


y - G (45) 

Substituting  (45)  into  (44),  the  following  result  is  ob- 
tained: 


S . . 
i J 


= G( 


2 v 


1 - 2 v Ykk 


6.  . 

i J 


2 y . . ) 
i J 


(46) 


Ykk  Yyy  + y3  3 


Also, 
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Equation  (46)  can  also  be  written  as: 

Sa3  = Gl-  T^Iv  (yyy  + Y33)  6aB  + 2YaB] 
Sa  3 2Gya3 


S 


33 


2G 

1 -2v 


+ ( 1 v ) Y33l 


(47) 


In  the  faces,  the  transverse  normal  stresses  are 
negligibly  small,  i.e.,  - 0 which  gives,  by  (47), 


and  when  this  is  substituted  in  the  first  set  of  equations 
(47),  the  following  result  is  obtained: 
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a3 
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(48) 


For  each  layer.  Equations  (17)  are  first  substi- 
tuted into  the  strain-displacement  relations  (3),  the 
Equations  (3)  are  then  substituted  into  (48)  to  get  the 
expressions  for  the  stresses,  and  finally  these  stresses 


are  substituted  into  (24)  to  get  the  stress-resultants 
in  terms  of  the  displacement  functions.  The  results  are: 
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(49) 

Substitution  of  Equation  (49)  into  Equations  (38) 
through  (43)  results  in  obtaining  a complete  formulation 
of  a general  and  nonlinear  problem  of  sandwich  plates  in 
terms  of  displacement  functions. 


CHAPTER  3 


STABILITY  ANALYSIS  OF  ELASTIC  SANDWICH  PLATES 

Combined  Axial  and  Transverse  Loads 
and  the  Problem  o f Buckling 

A special  case  of  the  analysis  presented  in  Chapter 
2 will  be  treated  here.  0A  is  now  considered  to  be  con- 
sisting wholly  of  0As  alone,  and  the  given  surface  forces 
over  0As  are  all  acting  perpendicularly  to  the  plate. 

Also,  the  body  forces,  moments  and  the  acceleration  forces 
and  moments  are  negligible  compared  to  other  forces.  Then, 


0 


and 


0 (50) 


The  axial  stress  resultants  N'  , N".  and  N _ 

ap  ap  aB 


can 


be  expressed  as: 


50 


N " „ = N°q  + n " „ 

aB  aB  aB 


J D = N D + n 0 

aB  aB  aB 


(51  ) 


where  N°  and  N D are  initial  axial  stresses  in  the  face 

aB  aB 

layer s and  in  the  core  respectively  before  transverse  de- 
flection occurs,  and  n'  , n''  and  n 0 are  additional  small 

aB  aB  aB 

stresses  in  the  upper  face,  lower  face  and  the  core  re- 
spectively, giving  rise  to  buckling. 


It  is 

assumed 

that : 

n 1 _ , 

aB 

« Na6 

and 

"o6  ""  fa6 

(52) 

i • e . , 

n ' and 

aB 

n1'  are 

aB 

much 

smal 1 er 

than  N 

° and  n _ are 

aB  aB 

much 

smaller 

than  N _ . 

aB 

The 

va 1 ue  s 

of  n ' 

aB 

and  n ''  are 

aB 

equal 

to  the 

values  of 

N ' 

aB 

and  N '' 

aB 

respectively  after 

drop 

ping 

all  the 

nonlinear 

terms . 

From 

( 50 ) and 

(38)  , 

the  following 

equilibrium 

equa 

ti  ons 

are  obtai ned : 

N 1 + — (M1 

a3,3  tf  ' aB , B 

- Q') 
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= 0 

(53) 

NaB,B  ' 

— — fM" 
tf  uaB,B 

- Q") 
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= 0 

(54) 

-r-  (Q1 

t^r  V^a 

- M 1 

aB,B 

- Q"  + 

a 

M"  ) + 

aB , B 

NaB,B 

0 

(55) 
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~ (Q*  - M'  a + Q"  - M"  J + M c Q - Q = 0 

a a3,3  a a3,B  a3>3  a 


(56) 


{ (N 1 q+N"q+N  JlL  q } + { (N  ' -N"J  -4  V.  . + M J 
a3  aB  a3  3,3  ?ot  a3  aB  2 3,3  a3  3,p,a 


<\j 


'Xj 


+(«;  + «;  + v,a  + (s3*>h+<s3*>-h'° 


(57) 


tr  t 2 tr 

4((N'  -N1'  ) U,  „}  + -4-  (N'+N"  )Y,  Q}  +-£(Q,-Q") 

2 a3  a3  3,3  ,a  4 aB  a3  3,3  ,a  2 a a ,a 


t ^ 

+ (M  qUq  q+K  'I'  0+T  ) - N„  + {(S'  ).  - (S"  ) , } = 0 

a3  3,3  a3  3,3  a ,a  33  2 3*  h 3*  -h 


(58) 


From  (53)  and  (54),  the  following  relations  are  obtained: 


-TT-  ( Q 1 - M'  „ ) = N'  Q 

t^  a aB  ,3  a3 , 3 


t---  (M"q  ~ Q"  ) = N"0 

t ^ a8 , 3 a aB  , 


(59) 


Then  (59)  is  substituted  into  (55)  and  the  result  is: 


N'  +N"  +N  =0 

aB,3  a3,3  aB,B 


(60) 


Substitution  of  (59)  into  (56)  gives: 


— I M _ n ) + ( N ' - N " ) 

tc  ' aB,3  u aB  aB  ,3 


0 


(61) 
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Using  Equation  (60)  in  (57)  the  following  equation  is 
obtained: 

t 

N'  + N1'  + N .)  U,  0 + { (N ' - N"J  4 + M 'i'  D} 

a6  aB  aB  3,Ba  aB  aB  2 3,B  aB  3,B  ,a 


+ % + + V,a  + <S3*>h  + %>-h  ’ 0 


(62) 


Subtraction  of  (54)  from  (53)  gives: 


t(M'  + M"  ) . - (Q'  + 0"  )]  + (N'  - N"  ) a = 0 


aB  aB  , B vva  va 


aB  aB  ,3 


(63) 


Addition  of  (53)  and  (54)  gives: 


[(M'  - M"  ) _ - (Q‘  - Q")]  + (N'  + N"  ) „ = 0 

t^  aB  aB  ,3  a a aB  aB  ,B 


(64) 


Now,  for  small  deformations  beyond  the  critical  point  in 
the  transverse  direction  the  middle  surface  of  the  core  is 
assumed  to  remain  undeformed, wh i ch  means: 


n a = 0 
aB 


Then,  from  (51)  it  follows  that 


NaB  = Nas 

Substituting  (51)  into  (61),  (63),  (64),  (62)  and  (58), 
and  making  use  of  (60),  the  following  equations  are  obtained 
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(M  Q Q - Q ) +(n'Q  - n"  ) Q = 0 

t aB  ,B  a v aB  aB  ,B 


(65) 


(66) 


t7[(Ma3  -MaB},B  - (Qd  ’ W + (r^  + n- ) jg  - ° (67) 


(Q'+Q"  +Q  ) + (2N°q+N  .)  U-  . + { (n ' -n"  ) -§-  Q + M Q) 

a a a ,a  aB  aB  3,aB  aB  aB  2 3,B  aB  3,3  ,a 


a.  a. 


+<S3*>h  + <S3*>  = 0 


(68) 


and 


t„2 


t„2 


2 {^naB  PaB^  U3,3^,a  + 4 ^2NaB  +naB+na3^ ,0^3,3  + T^a&^S.aP 


t ? 


+ 4 'naB+naB^  *3,aB  + 2 ^a"Qa\a  + ^aBU3,B+KaBV3,B+Ta^  ,a 


’N33  + 2 ^S3*^h  ^S3*^-h}  0 


(69) 


After  dropping  all  the  nonlinear  terms  from  Equations  (68) 
and  (69)  and  making  use  of  (60),  the  Equations  (68)  and 
(69)  become: 


«V.a+  <2I£b4V  U3,aB+  <S3*>h+  <S3*>-h  * 0 <70> 


and 
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)" ) + t N°  V 

a ,a  c a3  3,a3 


MS^)h  - (s$j_h  = 0 (71) 


In  a buckling  problem  in  which  deflections  in  the 


transverse  direction  are  small,  the  displacements  of  the 
middle  plane  of  the  core  in  the  direction  of  x can  be 


neglected  because  the  said  displacements  are  very  small 
compared  to  their  values  just  prior  to  the  occurrence  of 
buckling.  Therefore, 


This  reduces  the  number  of  unknown  quantities  from  ten  to 
eight  for  the  buckling  problem.  Consequently,  there  re- 
mains no  need  to  consider  all  the  ten  equations  given  by 
(55),  (65),  (66),  (67),  (70),  and  (71).  The  two  equations 
given  by  (55)  can  be  shown  to  be  redundant  and  are  there- 
fore dropped.  The  remaining  equations  are: 


a 


U 


0 


(72) 


a 


(73) 


(76) 
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and 


y=-  (T  - 
tc  a, a 


N33^ 


+ (Q'-Q")  „ + 


a a ,a 


2N°  W , Q 

a3  a3 


+ (p-  - £'•)  = 0 


where 


and 


w = u3 


^S3*^h 
P."  = (s"*)_h 


(77) 


(78) 


Here,  w is  the  transverse  deflection  of  the  middle 
surface  of  the  core  and  is  therefore  the  antisymmetric  or 
overall  deflection  of  the  sandwich  plate.  W is  the  trans- 
verse deflection  of  either  face  from  the  middle  surface 
of  the  core,  and  is  therefore  the  symmetric  deflection 
of  the  sandwich  plate.  The  quantities  £/  and  jd " are  the 
loads  acting  perpendicularly  to  the  upper  and  lower  surfaces 
respectively  and  are  positive  when  both  of  them  act  in  the 
x3  direction. 

The  simplified  Equations  (73)  through  (77)  were 
obtained  by  Kim  [42]  in  his  linear  analysis  for  the  com- 
bined axial  and  transverse  loads  problem.  The  present 
analysis  is  therefore  a general  treatment  of  which  Kim's 
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equations  can  be  obtained  as  a special  case. 


Kim's  Solution  of  the  Problem  of 
Combined  Axial  and  Transverse  Loads 


The  material  presented  from  here  on  up  to  the  end 


of  this  chapter  is  a summary  of  results  obtained  by  Kim 
[42]  in  deriving  the  expressions  for  the  antisymmetric  and 
symmetric  deflections  (w  and  W)  of  an  elastic  sandwich 
plate  from  the  Equations  (73)  through  (78).  These  expres- 
sions for  w and  W for  an  elastic  sandwich  plate  will  then 
be  used  to  derive  the  time  dependent  expressions  of  the 
same  two  quantities  for  a viscoelastic  sandwich  plate  making 
use  of  the  Correspondence  Principle  of  the  Viscoelastic 
Theory . 

The  following  quantities  are  defined  similar 
to  the  notations  used  in  [42]  and  [32]: 


u 


a 


t 


u 


ii 


a 


2 a 


t 


u 


a 


6 7 


f = - — N° 

ag  Gftf  aB 


q 2G 


V (£'  + £") 

ff 


2Gflf 


(p1  - p") 


Gctc 
Gf  tf 


and 


'Xj 

u 


a 


2 ( U 


a 


U" ) 

a 


/s 


(79) 


% 

First  u1  and  u"  are  solved  in  terms  of  u , u and 

a a a a 

the  results  are  substituted  into  the  expressions  for  V 

and  H'"  in  the  above  equations.  Then  the  various  stress 

a 

a> 

resultants  given  in  (49)  are  expressed  in  terms  of  the  u , 

A A A A 

u , u , f R,  f,  q,  q,  g,  w and  W defined  in  (78)  and  (79). 

The  non-linear  terms  and  the  quantities  uq  are  then  dropped 
from  these  expressions.  Keeping  in  mind  that  the  quantities 
n^  and  n^  are  given  by  the  linearized  expressions  of 
and  N ^ p respectively,  Equations  (49)  are  then  substituted 
into  Equations  (73)  through  (77).  Thus  the  equilibrium 
equations  in  terms  of  the  displacement  functions  and  initial 
stress  resul  tants  are  obtained  as  follows: 
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Hv 


6 < l^r  WVbb'  - t;  V%>  + Txf  <WW 


+Ua,BB  + Ua,BB  ^ 


(80) 


1+V.p  n 0/  4 ^ 9 2 ^ 

T^7  (uB,Ba+  3 W + ua,BB  + 3 Ua  ,BB  ' ^ t?  ““  + = ° (8'  ’ 


a. 


— — U + W,  + i-  (x^  U + W,  ) + f QW,  a + CJ  - 0 

t^  a ,a  aa  2 v tc  a ,a  aa  aB  aB 


% 


% 


(82) 


1-v^  UB,Ba  + Ua,B3  " 2tf  ^tf  Ua  + W’a^  0 


(83) 


(x~  U W , ) + r-  W , 

t^r  a, a aa  o 


4g(l-vc) 
a<X  tc2(1-2vc) 


W + f ,W  + q = 0 

LA  p Lip 


(84) 


There  are  eight  equations  and  eight  unknowns  in 
(80)  through  (84);  hence,  they  can  be  completely  solved 
with  the  appropriate  boundary  conditions.  For  a rectangu- 
lar plate  simply  supported  on  all  sides,  the 
solution  is  given  in  [42]. 

However,  in  a sandwich  plate  problem,  it  is  usually 
of  interest  to  calculate  w and  W only.  For  that  reason, 
eliminating  u and  u a from  Equations  (80),  (81)  and  (82), 
a single  equation  in  w is  obtained.  Similarly,  eliminating 
ua  from  (83)  and  (84),  a single  equation  in  W is  obtained. 
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They  are  as  foil ows  : 


R,V6w  + R9V4w  = (1-R3V2  + R4V4)‘  [(2N°a  + N„fl)  w,„R  + £'+£"]  (85) 


a8  a8  ’ cx3 


and 


4Df  (1  + f)  V4W  - G|-tc  [1  + . T ' C-: ] " + t n-'2V  V 

f 6 6 t 2 l-2vj(l-vj  V 4V 


i6ty(i-vc) 


2.,  , 4Gc'1-vc> 


C ' C' ' C' 


^1_  yrTv.T  ^ ^ ^2Na3W,a6  + ~ ^ 


(86) 


where 


2g(l-vc)(l-vf) 


DfD--  ^ 1 + f)  ^ + 3Tl-2v  ) 


G t 
c c 


g ( 1 - vc ) ( 1 - vf ) 


2tf  4 1 f 

D)[1  + .„_+  — Brrr^y-1-  ] 


D,  2tf  (l-vr)(l-vf) 

vr [1  + 9 ‘^7 + “ 6(,-2v  1 ] 


DfDi 


4 GfGctftc 


[1  + 


2g(l -v  ) ( 1 -vf ) 
3Tl-2vc) 


] 
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12(1 -vf 2 ) 


and 


= Young's  Modulus  of  the  face  layers. 


(87) 


The  order  of  the  differential  Equation  (85)  is  six, 


and  that  of  (86)  is  four,  therefore  these  two  equations 
require  three  conditions  and  two  conditions  respectively 
at  each  side.  Since  the  number  of  available  boundary  con- 
ditions for  w and  W is  one  each  at  each  of  the  four  sides 

of  the  panel,  in  general  there  are  not  enough  boundary  con- 
ditions to  use  Equations  (85)  and  (86)  independently. 


a and  b along  x and  y directions  respectively  is  considered. 
Let  x = x i , y = and  z = X3.  The  origin  of  the  coordinate 
axes  is  located  at  one  of  the  corners  of  the  middle  surface 
of  the  core.  The  plate  is  subjected  to  arbitrary  trans- 
verse loads  p'  and  p",  and  a uniformly  distributed  axial 
compressive  load  N along  the  face  edges  x = 0 and  x = a. 

Then  , 


Rectangular  Sandwich  Plate 
Simply  Supported  on  All  Si d e s 


A rectangular  sandwich  plate  having  sides  of  length 


N 


O 


N°  = N 

y x 


N 


0 


(88) 


xy 


y 
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where 


r = N , N°  = N , etc. 
x xx’  y yy  ’ 


Let 


00  GO 


w(x,  y)  = Z.._,  /L — i Sin  K,x  Sin  K0y 

m=l  n=l  mn  1 2 


CO  CO 

W Sin  K-.X  Sin  K0y 
m=l  n=l  mn  1 2 


CO  CO 

\ , T V - 

q(x,  y)  = Z . Z_i  o Sin  K,x  Sin  K9y 

m=l  n=l  Wn  1 Z 


and 


CO  CO 


i(x. 


m=l  n=l 


qm,_  Sin  K,x  Sin  K9y 
'mn  I 2 


(89) 


where 


v _ mir  i/  _ mir 
K1  a ’ K2  b 


mn  ab 


a b 

77  / q(x,  y)  Sin  K-j x Sin  K9y  dx  dy 


and 


a b 


'mn  ab 


ci ( x , y)  Sin  K^x  Sin  K^y  dx  dy 


(90) 
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If  Equations  (89)  are  substituted  into  Equations 


(85)  and  (86),  they  are  found  to  be  satisfied.  The  bound- 
ary conditions  for  various  cases  have  been  discussed  elab- 
orately in  [42]  and  for  the  present  case  it  has  been  shown 
[42]  that  if  the  expressions  given  in  (89)  for  w and  W are 
used,  all  the  boundary  conditions  are  taken  care  of. 


obtained  by  substituting  (89)  in  Equations  (85  and  (86), 
and  are  given  as: 


The  values  of  the  coefficients  w and  W 

mn 


mn 


are 


0 


[1  + B3  (Kj2  + K22)  + B4  (k/  + K22)  ] 


■mn 


(91) 


w 


mn 


mn 


and 


(92) 


mn 


where 


and 


73 


2 2 

aft/t  a 

1 = 6q  ^ 1 + 2^  ^ f + 9 “3^ 


6g 


^"r  4 2 9 ar 

B2  = T~  [af(1  + 2r  + | r ) + ~JJ  ] 


‘c2 


9 9 ga 

b3  ' ir  [“f  0 + 3 9r^  + -rf ] 


2 2 

aftr  tf  ga 

B4  ■ ^lg—  <“f  + -T-  > 


a f ^ f a 

B,  = (2  + I ) 


B,  = 2 (1  + 6a,  a r2) 
v f c ' 


6 6 


2 ga 


2 , 2 
3 aflf 


Gftf 


af  1 -v. 
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2(]-vc) 

“c  1 -2  vc  (95) 

The  above  results  will  be  used  in  the  next  chapter 
where  the  problem  of  a viscoelastic  sandwich  plate  will 
be  treated.  The  analyses  presented  at  the  end  of  the  pre- 
vious chapter  and  in  this  chapter  are  for  elastic  sand- 
wich plates  made  of  three  layers.  However,  the  procedure 
would  be  quite  similar  for  plates  made  of  more  layers. 
Also,  the  results  obtained  here  can  be  readily  reduced  to 

the  case  of  a plate  made  of  a single  layer,  thin  or  moder- 

% 

ately  thick.  If  all  the  primed  quantities  except  ( S ^ ^ 

% % 

and  ( S^* ) are  dropped  and  the  quantities  (S^^and  (S^*) 
are  replaced  by  a load  acting  in  the  transverse  direction, 
the  Equations  (73)  through  (77)  are  reduced  to  those  for 
a single  layer  plate  given  in  [54],  and,  in  addition,  if 
the  terms  ( S ^ * ) ^ and  ( S ^ * ) are  dropped,  the  equation  for 
buckling  of  a thin  plate,  as  given  in  [55],  is  obtained. 


CHAPTER  4 


BENDING  AND  STABILITY  OF 
VISCOELASTIC  SANDWICH  PLATES 


Method  of  Approach 

A direct  method  of  formulating  the  problem  of  any 
viscoelastic  body  (plates,  discs,  shells,  etc.)  is  to 
start  from  the  equations  of  motions  of  a continuous  medium, 
the  compatibility  conditions  and  the  constitutive  equations 
applicable  to  the  particular  viscoelastic  material  considered. 
Many  authors,  [46],  [43],  [47,  [56]  and  others,  have  used 
this  approach  for  plate  and  plate- type  problems. 

There  remains,  however,  another  method  to  obtain 
the  equations  or  solutions,  in  a rather  indirect  way,  of  a 
large  class  of  problems  in  linear  viscoelasticity  from  the 
known  equations  or  solutions  of  the  similar  problems  in 
linear  elasticity.  This  is  the  correspondence  principle 
which  will  be  explained  in  the  next  section.  This  method  of 
finding  the  solution  in  linear  viscoelasticity  is  especially 
useful  and  convenient  when  there  already  exists  a known  solu- 
tion to  the  "correspondi n g " elastic  problem  having  the  same 
geometric  and  boundary  conditions  (with  certain  restrictions). 
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Since  in  a classical  field  like  linear  elasticity,  thou- 
sands of  problems  have  been. solved  and  their  solutions 
are  readily  available,  they  can  be  used  to  get  the  solu- 
tions of  the  "corresponding"  viscoelastic  problems.  This 
later  method  has  been  used  by  [57],  [58],  [59]  and  others 
to  solve  plate  problems  and  various  other  types  of  prob- 
lems in  viscoelasticity.  In  the  present  analysis  also, 
this  latter  method  will  be  used  to  get  the  solution  of  a 
viscoelastic  sandwich  plate  problem  from  the  known  solu- 
tion of  the  "corresponding"  elastic  sandwich  plate  prob- 
lem which  has  been  given  in  the  last  chapter. 


Correspondence  Princi pie 


This  principle  has  been  stated,  derived  and  ex- 
plained by  Lee  [60].  However,  the  version  given  by  Pipkin 
in  [61]  of  the  same  principle  slightly  differs  from  that 
given  in  [60].  According  to  Pipkin,  the  one  given  in  his 
book  is  the  exact  version  of  the  correspondence  principle. 
Though  it  can  be  easily  verified  that  both  the  versions 
give  the  same  result,  the  one  given  in  [61]  will  be  stated 
below  and  used  thereafter  in  the  present  analysis. 


If  F = f(s)  = J e F(t)dt  = Laplace  transform  of  F 

O 

( t = time) 

and  F = s F , 


(96) 


"7  *7 
/ / 


it  has  been  shown  in  [61]  that  the  constitutive  equations 
and  the  equations  of  equilibrium  of  a linearly  viscoelas- 
tic material  take  the  following  form  in  the  transformed 
plane: 


a.  ■ = c 
i J 


i j kl  Uk  , 1 


and 


a . ■ ■ 

i J , J 


-pfi 


(97) 


where 


and 


a . ■ = stresses 
i J 


cijkl  ~ ^me  dependent  elastic  constants 


l).  = displacements 


density 


f ^ = body  force  components 


Now,  the  equations  of  an  elastic  medium  are 


a.  • - C . • i -i  U , i 

l j i J kl  k , 1 


and 


a.  . • 

i J 


p-f\- 


(98) 
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Comparing  (97)  and  (98)  it  is  found  that  the  s- 
multiplied  Laplace  transforms  of  the  time  dependent  visco- 
elastic quantities  satisfy  the  time  independent  elastic 
equations  for  a quasi-static  case  (where  inertia  effects 
are  neglected  due  to  slow  motions). 

Suppose  there  are  two  bodies  of  the  same  size 
and  with  the  same  boundary  conditions,  but  one  of  them  is 
elastic  and  the  other  one  viscoelastic,  and  for  a given 
set  of  loading  and  displacement  conditions  the  solution 
to  the  problem  of  the  elastic  body  is  known.  The  solution 
to  the  problem  of  the  viscoelastic  body,  which  is  now 
acted  by  a set  of  loading  having  the  same  spatial  distri- 
bution as  that  of  the  elastic  body  and  changing  with  time 
with  other  quantities  like  displacements,  elastic  moduli 
etc.,  is  obtained  as  follows: 

The  s-multiplied  transforms  of  all  variables  and 
elastic  constants  are  substituted  into  the  elastic  solu- 
tion and  the  inversions  are  performed. 

As  pointed  out  in  [62],  the  correspondence  prin- 
ciple is  not  applicable  in  the  following  cases:  (1) 

A point  on  the  boundary  of  a viscoelastic  body  may  be 
kept  free  of  stress  while  loads  are  being  applied  else- 
where and  after  some  time  a device  may  be  applied  which 
makes  this  point  undergo  a prescribed  displacement.  Then 
neither  the  surface  traction  nor  the  displacement  is  known 
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at  all  times  t>0,  and  neither  of  these  functions  can 
be  subjected  to  the  La  pi  ace  - tra ns f ormat i on ; (2)  The 

boundaries  of  a viscoelastic  body  may  be  changing  in  the 
course  of  time.  Melting  bodies,  bodies  solidifying  from 
melts  and  bodies  extruded  from  orifice  belong  to  this 
class. 


Deflection  of  Viscoelastic  Sandwich  Plate 
Under  Combined  Axial  and  Transverse  Loads 

In  order  to  use  the  results  (91)  and  (93)  for 
solving  the  corresponding  viscoelastic  problem,  the  first 
step  would  be  to  consider  all  the  time  dependent  quanti- 
ties and  separate  them  from  the  rest  so  that  finding  their 
Laplace  transforms  would  be  easy.  Apart  from  the  loads 
and  displacements,  other  time  dependent  quantities  are 
the  elastic  constants  v and  G.  However,  if  the  bulk  modu- 
lus K and  the  shear  modulus  G have  the  same  time-dependence, 
v is  constant  [61].  In  the  present  analysis  it  would  be 
assumed  that  v is  a constant  in  time  partly  because  this 
is  a fairly  valid  assumption  for  many  viscoelastic  mater- 
ials and  partly  in  order  to  make  the  analysis  simpler. 

The  assumption  that  visa  constant  has  been  made  by  many 
authors,  e.g.  [47],  [43],  [57]. 

Thus  in  Equations  (93)  through  (95)  the  two  quan- 
tities g and  f turn  out  to  be  time  dependent,  and  therefore 
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the  expressions  for  B-j  through  Bg  given  in  (95)  are  re- 
written as  follows  after  some  rearrangements  of  terms: 


and 


B1  g + A29  + A3 


A4  + A59 


- + A7 

g 7 


— + A 

g a9 


A1 0 + A1  1 9 


B6  = A1 2 9 


A1  3 9 


B8  A1 5 


where 


2 , 2.  2 
af  tf  tc 


'1  8 


A, 


2 2 

axa  t tjr 

f c c f 
18 


(99) 


2 2 

afactf  tc 


2afV 
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2 2 

t a ar 

_L_f L_  (_c  + _±\ 

c 1 3 2 ' 


aftf 


1 1 


A, 


itc2af(l  + 2r  + | r2) 


1 + 8a , a r‘ 
f c 


'12 


A5  lc  ac 


2a 


A 


1 3 


_c 

2 


A6  lc  af 


A, 


14 


2 , 2 
3 aftf 


? ? ? t 2a 

A7  3 tc  r af  12 


Let  the  transverse  loads  have  the  following 


R*  = Pi (x,y)F(t) 


and 


£"  = P2(x,y)F(t) 


where 


and 


F (t)  = 0 for  t < 0 


F (t)  = an  arbitrary  function  of  time 


Then  from  ( 79 ) , 


q = Hi  (~)  and  q = H^^) 


(100) 

form : 


(101) 


(102) 


(103) 
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where 


Pi  - P2 


and  H, 


2 1 


P1  - p2 
2U 


0 04) 


Substitution  of  (103)  into  (90)  gives: 


q = (tO1) 

^mn  3 'G^ 


and 


Pmn  H4 


(105) 


where 


a b 

4 f f 

H0  = r~  I H,  Sin  K,x  Sin  K,,y  dxdy 
3 ao  I 1 1 

J 


and 


a b 


4 a b 


Sin  K-jX  Sin  K^y  dxdy 


005) 


Now  the  Equations  (99)  are  first  substituted  into 
(93)  and  (94),  then  (93),  (94)  and  (105)  are  substituted 
into  (91)  and  (92)  to  get  the  expressions  for  w^^  as 
f ol 1 ows  : 


w 


mn 


h3f  (H5g  + V 

Gf{H7  + H8g  + Hgg2  - (H1Qg  + H^)  f} 


(107) 
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and 


, . ' hijlF ' 

mn  ^]4  ~ H i g 9 ) 


(108) 


where 


H5  1 + A7K3  + A9K3 


H6  A6K3  + A8K3 


H7  A1K3 


H = A K + A K 
8 3 3 


'43 


H 9 A2K3  + A5K3 


H10  = (1  + A7K3  + A9K32)  Kl2 


Hn  = (a6k3  + a8k32)  K]2 


^ ] 2 ^4  ^ + A14^3^ 


H 1 3 A11K3  + A 1 2 K 3 + A13 


H 1 4 A1 0K3 


H 1 5 ^ + A 1 4 K 3 ^ K1 


and 


K3  K1  + K2 


009) 
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After  substituting  g = anc'  f = 

into  Equations  (107)  and  (108),  w and  W are  obtained 
^ v ' ' mn  mn 

in  the  following  form: 


w 


FH 1 6 (Gc  + H17Gf} 


{Gc  + H18GcGf  + hl  1 9 G f “ ^ H20Gc  + H 2 1 G f ^ 


(no) 


and 


F H 


W 


37 


mn 


Gc  + H 3 8 G f 


h39n 


(111) 


w h e r ! 


H 0 H r-  r 

O D 


'16  H 


H?r 

I - _L_ 

19  H n 


H r 

I = - ■ ■■■■■ 
1 1 7 H r 


^1  or 

n = Lr 

'•20  Hgtf  ’ 


H18  - -f 

H9 


Hnr< 


21  Hgtf 


‘37  H 


U 
1 3 


38  H 


11 

1 3 


and  H 


39  H 


15 

1 3 


(112) 


The  Equations  (110)  and  (111)  are  for  an  elastic 
sandwich  plate.  According  to  the  correspondence  principle 
discussed  in  the  last  section,  these  two  equations  for  a 
viscoelastic  sandwich  plate  would  take  the  form: 


(113) 


w 

mn 


H1 6 <Gc  + H1  7 Gf ^ 


5c?  + H,8GcGf  + HT9Gf2  - <H20Gc  + 


H2lGf)N 


and 


W 

mn 


F H 


37 


H38Gf 


(114) 


At  first,  the  deflections  will  be  found  for  F(t) 
equal  to  a unit  step  function  and  it  will  be  shown  that 
the  solution  thus  obtained  could  be  used  for  any  arbitrar- 
ily time  dependent  load  with  the  help  of  convolution  in- 
tegrals. Then,  for 

F(t)  = H(t),  the  Heaviside  unit  step  function, 

F = ? 015) 


Also,  the  axial  load  is  assumed  to  be: 


N ( t ) = N.  H ( t ) 

¥ = ^ i (ii6) 

s 

For  the  values  of  F and  N given  in  (115)  and  (116). 
Equations  (113)  and  (114)  now  become: 
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w 


mn 


NT  <"7) 


and 


W 


mn 


(118) 


Materia]  Parameters 

It  is  assumed  that  the  core  is  made  of  an  isotropic 
and  linearly  viscoelastic  material  w rich  can  be  represented 
as  a 3-parameter  model  as  shown  in  Figure  (4).  The  upper 
and  lower  face  layers  are  assumed  to  be  made  of  the  same 
materia],  but  different  from  that  of  the  core.  The  mater- 
ial which  the  faces  are  composed  of  is  also  assumed  to  be 
isotropic  and  linearly  viscoelastic  represented  by  another 
3-parameter  model  also  shown  in  Figure  (4).  Each  of  the 
two  models  has  a spring  in  series  with  a parallel  combina- 
tion of  a spring  and  a dash  pot  and  therefore  represents 
a viscoelastic  solid.  From  here  on,  the  subscripts  c and 
f will  refer  to  the  core  and  the  face  materials  respectively. 

From  the  tabulated  chart  given  in  [62],  the  value 
of  the  relaxation  modulus  in  shear  for  the  core  material 
is  given  by: 


87 


y 


c2 


ycl 

AW- 


• o-AW — 


( core ) 


Figure  4 


■°-W\/V — 


y 


f 2 


yfl 

-4A/W 


f 

( faces ) 

Th ree- parameter  models  for  viscoelastic 
materials 


Y (t)  = -p--1  e 
c cl 


p , p 

^el  + pc2  " e c ^ ' 


(119) 


where 


and 


ycl  + yc2 


c2 


y c ! + y c 2 


c3 


ncyc2 


ycl  + Uc2 


(120) 


From  (119)  and  (120)  it  follows  that  for  the  core. 


Pc6<s  * Pc5> 
5(5  + Pc4) 


(121) 


oo 


where 


c4 


c 1 


c5 


c 2 
c 3 


P = P P 
c6  c3  c4 


In  a similar  way,  the  expression  for  the  relaxation  modu- 
lus of  the  face  material  can  be  obtained  as  follows: 


Yf(t) 


>3 

fl 


e P 


fl  + P f 2 0 - e Pfi  ) 


(122) 


and  therefore. 


= Pf6  (S  + Pf5) 
s(s  + Pf 4 ) 


(123) 


in  which  the  quantities  Pf]  through  P^-g  are  defined  in  a 
way  exactly  similar  to  that  for  the  core. 

Since  in  the  present  analysis  the  external  loads 
applied  to  the  plate  are  Heaviside  unit  step  functions  of 
time,  the  quantities  Gc  and  in  Equations  (117)  and  (118) 
are  the  same  as  the  quantities  Yc  and  Y^  respectively  in 
Equations  (119)  and  (122). 

Antisymmetric  Deflection 

In  Equation  (117),  the  quantities  Yc  and  Y ^ whose 
values  are  given  in  (119)  and  (123)  are  substituted  for 
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Gc  and  respectively,  and  the  result  is: 


w 


mn 


H 30  ^ s + H31S  + H32 ) ( S + PC4 ^ s + p f 4 ) 
s(s  + H33sv  + H34s  + H35s  + H36) 


024) 


where 


H 1 6 H 3 2 


'2  3 


'30  H 


25 


31  H 


22 


24 

32  H22  ’ 


26 


'33  H 


25 


'34 


27 
H2  5 


'28 


'29 


'35  H 


25 


'36  H 


25 


'22 


(Pc6  + H17Pf6> 


H23  ~ Pc6^Pc5  + Pf4)  + H17  Pf6(Pf5  + Pc4) 


H24  Pc6Pc5Pf4  + H17Pf6Pf5Pc4 


H25  Pc62  + H18Pc6Pf6  + H 1 9 P f 6 ' ^H20Pc6  + hl  2 1 P f 6 ^ N i 


H26  ' 2Pc6  (PC5  + Pf4>  + H18Pc6Pf6(Pc4+Pf4+Pc5+Pf5>  + 2H19Pf6<Pc4+Pf5) 


-(H2uPc6(2Pf4+Pc5+Pc4>  + H21Pf6<2Pc4+Pf5+Pf4>}  Ni 


H27  = Pc6<Pc5+p2f4+4Pc5Pf4>  +H18Pc6Pf6(Pc5Pf5+Pc4Pf4+<Pc5+Pf5)(Pc4+Pf4>} 

( con  td  . ) 
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,2  /n2  , n2 


+H,9Pf6<Pf5+Pc4+4Pf5Pc4>  ' [ Vc6{Pc5Pc4+tf4(Pc5+Pc4>  } 


+H21 Pf6  ! Pf5Pf4  + Pc4  + 2Pc4  *Pf5  + Pf4*}^  Ni 


»28  = 2Pc6<Pc5p2f4tPf4Pc5>  + H18Pc6Pf6{<PcS+Pf5)Pc4Pf4+(Pc4+Pf4>Pc5Pf5} 


+2H19p2f6(Pf5Pc4+Pc4P?5)  " [H20Pc6(<Pc5+Pc4>P?4+2Pf4Pc5Pc41 


+H21Pf6  {^Pf5+Pf4^Pc4  + 2Pc4Pf5Pf4  Ni 


H9o  = pLp2cPL  + H, 0P  CP,CP  ,P,CP  „P.,  + H,,,P2  P2  P2 
29  c6  cb  f4  18  c6  f6  co  f5  c4  t4  ;9  f6  fb  c4 


.(H  P P P P2  + H P p P P2  ) N 
1 20  c6  c5  c4  f4  M21  f6'  f5‘  f4  c4'  i 


(125) 


The  next  task  is  to  find  the  inverse  transform  of 
the  expression  in  (124).  As  the  first  step  to  achieve 
this,  the  right  hand  side  of  (124)  is  written  as  follows 
by  the  method  of  partial  fractions: 


w 


mn 


B i B^s2  + B^s2  + B^s  + Bp- 

S + H58S  3 + P 5 9 S ^ + H60S  + H61 


= H 


30 


r h.  + Q(s)  -i 

^ S P ( S ) 


(126) 


y i 


where 


Q(s)  = s ^ + B^s2  + B^s  + B5 


(127) 


and  the  quantities  B -j  , B 2 , B 3 , B ^ and  B ^ can  be  easily 
found  using  the  method  of  partial  fractions. 

Let  b-j  , b 2 > b3  and  b^  be  the  roots  of  the  poly- 
nomial P ( s ) , i . e . , 


The  inverse  transform  of  w r ^ can  have  different  ex- 


pressions depending  on  the  values  of  the  roots  b ^ , b 2 , b 3 
and  b . All  these  different  possibilities  will  now  be 
explored. 

Case  ( I ) : All  Real  Roots 

Within  this  category,  there  can  be  several  possi- 
bilities as  follows: 


From  the  table  of  Laplace  transforms  [63]  the  inverse 

transform  is  obtained  and  the  value  of  w is  given  by 

mn  3 J 


P(s)  = (s  - b-jMs  - b ^ ) ( s - b 3 ) ( s - b4) 


(128) 


( a ) All  roots  different  from  each  other 


i . e . , b -j  f b2  1 b3  f b4 
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where 

Ltli 

pk(s)  = s'bk 

( b ) Tw o equal  roots 
Let  bi  = b2  f b3  + b4 

Let 

Q ( s ) B6S  + B7  B8s  * B9 

PITT  ' (S  . b))2  (s-b3) ts-b4) 

in  which  , B 7 , Bg  and  Bg  can  be  found  easily. 

tuting  the  above  quantity  into  (126)  and  ref err 

the  following  result  is  obtained: 


wnin^^  H30 ^ B 1 + B6  {1  + ^B10  + bl^t}  e 


V 


bg+Bii  b q t b4.  + bii  bab 
+ B8  {-^—^1  e 3 + 4 11  e 4 }] 


b 3 b4 


b4'b3 


where 


and  B-j  -j 


( c ) Three  equal  roots 
Let  b]  = b2  - bg  f b4 


(129) 


Substi - 

inq  to  [26] 


(130) 


Then  from  (126) 


b 3 + B10bfl2  + B, ^b . + B- 


wmn^^  H30  ^ B 1 + B2  { 


(b4  " b]) 


b „ t 


r 4 1 2 4 1 3 4 1 4 4 


3b4b,2  - 3b42b,  - b 3 - B12b  2 - B)3  + Bh  b t 

e 

(b4  - b,)3 


( B1  2-3b4  ) bi  + 2bl  2 B 1 2 b4  b 1 ‘ B 1 3 b 4 B14  x b 

+ : — — - — 5 te 

(b]  - b )2 


B 1 2 b 1 2 + b]3  + B13bi  + B 1 4 9 b 1 1 3 

+ t e ) J 

2(b,  - b4) 


where 


B 3 B4  B 5 

D - D - 2.  D J. 

12  B2  ’ °1 3 B2  ’ D 1 4 B2 


(131) 


( d ) A1 1 equal  roots 

Here  b-j  = b2  = b^  = b^. 

Then  the  inverse  of  (126)  is  given  by 


B,  T + 2B-,  „b,  + 3b, 

wmn  (t)  ■ H3cfBl  + B2  <’  + <B12  + 3bl>t  + 1 


B14  + B 1 3 b 1 + B 1 2 b 1 + bl  3 b i b 

+ lAJ _ !—  t3}  e 1 ] (132) 


where  the  values  of  B ^ 2 , B^2  and  B^  are  given  in  (131). 
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( e ) Two  pairs  of  equal  roots 

In  this  case,  b -j  = bp  f b ^ = b 4 

This  can  be  treated  as  a special  case  of  all  complex  roots, 
to  be  treated  next,  in  which  the  complex  parts  of  the  roots 
are  zeros . 


Case  (II):  All  Complex  Roots 

Since  the  complex  roots  occur  in  pairs,  let 

b •]  = a -j  + i 3 1 

bp  = a -j  - ’ B-, 

^3  = a2  + ^2 

t>4  ~ ct  p " ^ B p 

Then  , 

(s  - b i ) ( s - bp)  = (s  - a-j  ) 2 + 3 -| 2 

and 

(s  - b 3 ) ( s - b4)  = (s  - a2)2  + Bp2 

Hence  , 

P(s)  = {( s - a-j  ) 2 + B -j  2 j {(  s - a2)2  + 3p2J 


Now,  let 


Q ( s ) 

PTT7 


D-j  s + D2 


D3s  + D4 


(s-a1 )2  + 31 2 


(s-a2)2  + Bp2 


+ 


(133) 


55 


in  which  the  quantities  D-j  , Dp , Dp  and  can  be  solved 
by  the  method  of  partial  fractions. 

Substituting  (133)  into  (126)  and  looking  up 
standard  tables  of  Laplace  transforms  like  [63], 
the  result  obtained  is: 


cut  D,cu+D„  out  a?t 

Wmn  ^ = H30  tBl  + Dle  cos  ^ + — g ~ e sin  3-jt  + D3e  L cos  Bpt 


D?a?  + D,  a„t 
+ g — e c sin  Bpt] 


(134) 


For  the  case  of  two  pairs  of  real  and  equal  roots 
as  mentioned  earlier,  the  results  given  in  (134)  car,  be 
used  by  taking  3-j  = Bp  = and  to  avoid  division  by  zero 

during  computation,  the  quantities  S'irK[lt  anfj  Sin c. ^ t can 

3 1 Bo 


S i n B -I  t SinB„t 

be  written  as  ( —D~ — )t  and  ( — — — ) t respectively,  and 


B-j  t 


Bp  t 


these  two  quantities  then  become  equal  to  t as  B^O  and 


Bp^-0 . 

In  a similar  way,  the  results  given  in  (134)  can 
also  be  used  in  the  case  where  there  are  a pair  of  complex 
roots  and  a pair  of  real  and  equal  roots. 


Case  III:  Two  Real  and  Unequal  Roots  and  a Pair  of  Complex 

Roots 


Let  b-|  = y + 


96 


and 


and 


Then  , 


b ^ f b ^ where  b ^ and  b ^ are  real  roots, 
(s  - b]  ) (s  - b2)  = (s  - y)2  + 62 . 


Let 

nr  \ Dcs  + Dc  D-.S  + D0 

Q ( s ) _ _5_ 6 7 8 

P ^ (s-y)2+S2  (s-b3)(s-b4) 


(135) 


where  the  quantities  D g , D g , D-,  and  D0  can  be  readily  found 
by  using  the  method  of  partial  fractions. 

Substituting  in  (126)  and  upon  performing  the  in- 
verse transform,  the  following  result  is  obtained: 


w 

mn 


(t) 


. D,-Y  + Dfi  . 

Hon  [B,  + DceYC  Cos  6t  + — — j e^ w Sin  6t 

3 U 1 b 6 


b + Dg  b3t  b4  + °9  b4tn 
+ °7(  F7-nT  6 + b7^-b7  6 )] 


(136) 


where 


Case  (IV):  Two  Real  and  Eg ual  Roots  and  a Pair  of  Complex 

Roots 

This  is  a special  case  of  (II)  where  either  8 -|  or 


$2  Is  zero . 


After  w (t)  has  been  found  in  any  of  the  above 
mentioned  ways,  it  is  to  be  substituted  into  the  series 
expression  given  in  (89)  to  compute  the  antisymmetric  de- 
flection w(x,  y,  t)  for  any  value  of  x,  y and  t. 

Symmetric  Deflection 


replaced  by  Yc  and  Y^  respectively  and  the  values  of  Yc 
and  Y ^ as  given  in  (121)  and  (123)  are  substituted.  The 
following  result  is  obtained: 


In  Equation  (118),  the  quantities  Gc  and  Gf  are 


W 


(137) 


where 


and 


(138) 
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Let  (s2  + H4gs  + H47)  = (s  + a)  (s  + $). 

The  following  three  cases  may  arise  depending  on 
of  a and  3 : 


Case  (I):  Real  and  Unequal  Ro ots 

Here,a7  8;aand  3 are  real. 

Then  taking  the  inverse  transform  of  (137),  it  is 
that 


W (t;)=  H45[  ^ e-at  _ B,  ~ H43B  + H44 

45  H47  a(a  - 6)  3(a  - 8) 


Case  (II):  Real  and  equal  Roo t s 


This  can  be  treated  as  a special  case  of 
be  considered  next. 


C a s e (_I JJ  ) : Comp  1 ex  Roots 


Let  a = y.|  + i6^ 


a = 


Y1 


i 6 


1 


Then  the  inverse  transform  of  (137)  is: 


W (t) 
mn  v ' 


H44  n "V  , D11  D10  Y1  'Y1 

‘45 L H„  D10  Cos  6lt  6,  e 


c;  C 


'47 


1 


where 


H 

H 


44 

47 


the  values 


found 

e"St]  (139) 
(III),  to 


t 

Sin  6 -j t ] 
(140) 


1 


99 


and 


1 1 


43 


H44H46 


'47 


The  result  (140)  can  be  used  for  the  case  where 
the  roots  are  real  and  equal  by  putting  6 ^ = 0.  Division 
by  zero  can  be  avoided  as  shown  in  connection  with  the 
antisymmetric  deflection. 

The  value  of  W (t)  so  obtained  can  now  be  substi- 
tuted into  the  expression  given  in  (89)  to  obtain  the  sym- 
metric deflection  W(x,  y,  t)  for  any  x,  y and  t. 

An  important  remark  can  be  made  at  this  point.  In 

whichever  of  the  above  mentioned  ways  w (t)  or  W (t) 

mn  mn 

are  calculated,  it  is  always  found  to  be  consisting  of  two 
parts,  a constant  term  plus  a few  exponential  terms.  For 
N . less  than  its  critical  value  for  either  mode,  antisym- 
metric or  symmetric,  for  any  given  t (it  will  be  shown  in 
the  next  chapter  that  the  critical  value  of  the  axial  stress 
resultant  N is  a function  of  time),  both  w (t)  and  W (t) 
must  be  finite,  and  this  requires  that  the  real  parts  of 
the  roots  of  P(s)  and  (s  + H ^ g s + H ^ ^ ) be  negative.  It 
is  easy  to  prove  that  this  requirement  is  satisfied  for 
(s  + H^gS  + H ^ ^ ) . It  is  seen  that  in  P ( s ) , all  the  co- 
efficients are  of  the  same  sign  which  is  a necessary  but 
not  a sufficient  condition  for  the  real  parts  of  all  the 
roots  to  be  negative.  To  see  that  both  the  necessary  and 


TOO 


sufficient  conditions  are  satisfied,  one  should  check  the 
Routh -Hurwi.tz  criterion  which  is  given  in  any  standard 
text  on  the  stability  of  dynamic  systems  like  [64].  How- 
ever, it  can  be  shown  numerically  (numerical  results  are 
given  in  the  last  chapter)  that  this  condition  is  indeed 
satisfied,  i.e.,  the  time  dependent  terms  in  the  expres- 
sions of  wmn(t)  and  W (t)  are  exponentially  decreasing 
functions  of  time.  A few  more  remarks  have  been  made  in 
the  discussion  chapter. 

Thus  , 


w 


mn  (at  t = °°) 


H 3 0 B 1 


and 


W 


m n (at  t = °° ) 


H45H44 


'47 


041) 


Stability  Analysis 

If  the  transverse  loads  p-j  and  p^  are  zero,  it  is 
found  from  (104),  (106),  (109)  and  (112)  that  the  quanti- 
ties H -j  g and  H ^ 7 are  also  equal  to  zero.  Then  from  (113) 

and  (114)  it  follows  that  the  deflections  w and  W of 

mn  mn 

the  elastic  sandwich  plate  can  grow  indefinitely  if  N takes 
up  the  following  two  values  correspondi ng  to  the  antisym- 
metric instability  and  the  symmetric  instability  respectively 
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cr 


(a) 


Gc2  + H1SGcGf  + H1 9Gf 2 
H20Gc  + H21 Gf  ' 


(142) 


(anti symme tr i c ) 


and 


cr  ( s ) 


Gc  + H38bf 
H 39 


(143) 


( symmetr i c ) 


where  the  quantities  Ncr(a)  and  ^cr(s)  denote  the  critical 
values  of  N for  antisymmetric  and  symmetric  buckling  re- 
spectively. 

According  to  the  method  used  earlier,  the  time 
dependent  critical  values  of  N for  a viscoelastic  sand- 
wich plate  are  obtained  from  (142)  and  (143)  as  follows: 


cr  ( a ) 


Sc2  4.  H18GcSf  * H19Gf2 
H20Gc  + H 2 1 Gf 


or , 


Ncr(a) 

and 


Gc  + H18GcGf  + H18Gf 
H20Gc  + H21Gf 


(144) 


N 


cr  ( s ) 


Gc  + H38Gf 


'39 


or 


cr  ( s ) 


G + H-nGf 
c 38  t 


39 


(145) 
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Replacing  G-  and  Gf  by  Yc  and  Yf  respectively  in  Equations 
(144)  and  (145),  and  substituting  (121)  and  (123)  for  Yc 
and  Tf,  the  following  results  are  obtained: 


H56^4  + H57S  + H58S  + H 5 9 S + H60j| 


cr(a)  ' s(5  = Pc4)(s  + Pfa)(s2  + HS1S  + Hfil) 


f 4 


'61 


‘61 


(146) 


and 

- = H66^S  + H67S  + H68  ^ 

Cr^  s(s  + Pc4) (s  + Pf4) 


(147) 


where 


'56 


"48 

^53 


'57 


49 

!48 


5 8 


50 

'48 


‘59  H 


5.1 

48 


u - 52  u = _54  u 

60  "H48  ’ 61  H5  3 ’ 62 

The  quantities  H48  through  H52  are  equal  to  the  expressions 
for  the  quantities  H25  through  H2g  respectively  given  in 
Equations  (125)  if  N . is  set  equal  to  zero. 


'55 


'53 


H53  = H20Pc6  + H2 1 Pf 6 


H54  " H20Pc6(Pc5  + Pf4)  + H2 1 Pf 6 ( Pf 5 + Pc4} 


H55  = H20Pc6Pc5Pf4  + H46Pf6Pf5Pc4 


'66  H 


6 3 
39 


'67  H 


64 

53 


H 


H 

68  = H 


65 

63 


H6 3 ' Pc6  + H38Pf6 


H64  " Pc6(Pc5  + Pf4)  + H38Pf6(Pf5  + Pc4) 


and 


H65  = Pc6Pc5Pf4  + H38Pf6Pf5Pc4 


0 48) 


The  inverse  transforms  of  (146)  and  (147)  have  to 
be  found,  as  will  be  shown  below. 

Antisymmetric  Buckling 

Equation  (146)  is  written  as  follows: 


cr ( a ) H56^  s + (s+P-4)(s+Pf4)  $2  + H s + h 


Bs  + C 


Os  -+  E 


- ] 049) 


c4 


61 


'62 


in  which  the  coefficients  A,  B,  C,  D and  E can  be  found 
easily.  The  inverse  transforms  of  the  three  terms  within 


the  brackets  will  be  found  separately. 

A 

(I)  Inverse  transform  of  j = A 

B s + C 

(II)  Inverse  of  transform  of  p ^ ) ( s + P~T 


(150) 


B + ^ 


( B ^ ” P f 4 


c4  " f 4 


P f 4 “ Pc4 


] • 


• if  Pc4  * P f 4 


= B [1  + (fe- 


B c4 


P.,)L]  e 


Pc4l 


if  Pc4  = P f 4 


(151) 


(Ill) 
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Let  s2  + H 6 1 s + H 6 ^ = (s  + <i>-,)(s  + <f>  2 ) 


(a)  If  and  4>2  are  real,  finding  the 
inverse  of  the  third  term  within  the  brackets  in 
(149)  becomes  similar  to  that  for  the  second  term 
in  the  same  expression. 

(b)  If  and  <t>2  are  complex,  let 
= 0 1 + i02  , and 

4>  o ~ - I ® 2 


Substituting  the  above  expressions  for  the  inverse 
transforms  in  (149),  the  value  of  ^cr(a)(m’  n>  t ) is  now 
obtained. 

S.ymmetr  i c Buckling 

From  Equation  (147),  the  following  results  are 
obtained: 


Then  the  inverse  transform  of  the  third  term 


D [e 


- 0 1 1 


Cos  0ot  + 


Sin  9 ^ 1 3 (152) 


f 4 6 7 f 4 68 

Pf4^  Pc4  " P f 4 ^ 


e 
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and 


H 


Ncr(s)(m’"’t)  = H66[, 


68 


H 


c4 


68_ 

>2 

c4 


i _Pr4t 

1 + (P,4  - Hfi7  + P-  > 0 6 ] 


c4 


for  P 


f 4 


(153) 


It  can  be  easily  shown  that  the  time  dependent 
parts  in  the  expression  of  both  Ncr(a)  and  Ncr(s)  are  ex- 
ponentially decreasing  functions  of  time.  Then  at  t = , 


and 


Ncr(a)(at  1 ' ”>  = H56  A 


„ , . . , H66H6S 

f!cr(s)(at  ‘ Pc4Tp: 


(154) 


Numerical  results  have  been  given  in  the  last  chap- 
ter, where  the  results  obtained  in  this  section  have  been 
discussed. 


Elastic  Faces  and  Viscoelastic  Core 

The  analysis  of  a sandwich  plate  with  elastic  faces 
and  viscoelastic  core  comes  out  as  a special  case  of  the 
analysis  presented  in  this  chapter.  In  the  three-parameter 
viscoelastic  model  of  the  face  layers  shown  in  Figure  (4), 
if  either  yf1  or  nf  is  set  equal  to  infinity,  or  if  nf 
if  set  equal  to  zero,  the  model  is  reduced  to  a single 
spring  of  stiffness  equal  to  u^2  or  two  springs  of  stiff- 
ness and  yf2  in  series,  either  of  which  represents  an 
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elastic  material.  However,  it  is  found  to  be  convenient 
if  nf  is  set  equal  to  infinity  which  means  that  P^  = °° 
and  = °°,  and  therefore  Pf4  = 0 and  Pf5  = 0.  But 

setting  P^  = 0 and  P^^  = 0 in  the  final  expressions  of 
wmn,  Wmn,  Ncr^a^  and  Ncr^sj  would  end  up  with  expressions 
like  zero  divided  by  zero.  To  avoid  this  difficulty,  one 
should  set  Pf4  = 0 and  Pf5  = 0 in  (124),  (137),  (146)  and 
(147)  and  calculate  the  quantities  w , W , N , , and 
Ncr(s)  respectively  following  a procedure  exactly  similar 
to  the  one  shown  before. 

Response  to  a L pad  wh 1 c h Is 
an  Arbi trary  Function  of  Time 

In  the  analysis  of  viscoelastic  sandwich  plate  pre- 
sented so  far,  it  has  been  assumed  that  the  transverse 
loads  p‘  and  p"  are  step  functions  in  time,  and  the  values 
of  wmn(t)  anc*  ^mn(t)  stained  are  therefore  corresponding 
to  step  loads.  But  these  results  can  be  generalized  for 
loads  which  are  arbitrary  functions  of  time  as  follows: 

Let  the  transverse  loads  p1  and  p"  be  uniform  through- 
out the  plate.  Then,  from  (106), 


1 6H, 
H3  - 2 ' 

1 6H  9 

m 

=1,3,5, 

and  = 

2 

=1,3,5, 

it  mn 

ir  mn 

n 

Thus , for 

unit  loads  , 

, Hg , 

and 

become 
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H 3 ( n ) ~ a n d H 4 ( | \)  ~ ~ p — . ( U ) signifies  unit  load 

7T  mn  v ' -rr  mn 

and  let  the  corresponding  deflections  which  are  propor- 
tional  to  H3  and  H4  become  wmn(u)  and  Wmn(u).  Let 

£'  = P-j  (x,y)f  (t) 

and 

P"  = P2(x,y)f (t) 

where 

f ( t ) = 0 for  t £ 0 

= an  arbitrary  function  of  time  for  t>0 

Then  it  follows  that  the  deflections  w , , and  W / , 

m n ( a r ) mnlar) 

for  loads  which  are  arbitrary  functions  [f(t)]  of  time 
are  given  in  the  form  of  convolution  integrals  as  follows: 

fl  ,f 

wmn  (ar)  = j wn,n(U)^  " t)  H?  dT  (155) 

O 

and 

t 

Wmn(ar)<t>  = _/W mnfU)^  ’ T>  37  dT  (,56> 

o 

where  x = time  variable  and  the  subscript  (ar)  denotes  ar- 
bitrary function  of  time 

An  important  remark  can  be  made  at  this  point. 
Unlike  the  transverse  loads,  the  axial  load  N appears 


non! i nearly  in  the  expressions  of  w (t)  and  W (t); 
therefore  generalization  similar  to  the  above  cannot  be 
made  for  N. 

Special  Cases 

The  case  of  elastic  faces  and  viscoelastic  core 
has  been  discussed  already.  The  deflections  of  a sand- 
wich plate  under  transverse  loads  only  can  be  readily  ob- 
tained from  the  expressions  of  w(t)  and  W(t)  by  setting 
N = 0.  Simpler  models  (like  the  two -pa rameter  Maxwell 
or  Kelv in-Voigt  models)  for  the  viscoelastic  materials 
can  be  treated  as  special  cases.  Finally,  the  elastic 
sandwich  plate  problem  itself  can  be  obtained  as  a special 
case  of  the  present  analysis  by  modifying  the  values  of 

ycl’  yc2  ’ V ufl’  yf2  and  nf‘ 

Comparison  with  Other  Theories 

As  mentioned  in  the  first  chapter,  very  few  reports 
on  the  study  of  viscoelastic  sandwich  plates  are  available. 
Malmstrong's  work  [44]  considers  time  dependent  creep  de- 
formation of  circular  sandwich  plates,  but  does  not  con- 
sider instability.  Also  the  deformation  is  due  to  creep 
beyond  the  yield  point  and  in  the  plastic  range,  and  there- 
fore is  quite  different  from  the  present  analysis  except 
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that  both  analyses  have  considered  time  dependent  quasi- 
static deflections. 

The  work  by  Chang  and  Fang  [43]  is  closer  to  the 
present  analysis.  They  have  considered  the  dynamic  ef- 
fects fora  transverse  step  load  and  then  for  an  arbitrary 
transverse  load.  However,  they  have  neglected  the  viso- 
elastic  effects  of  the  faces,  the  compressibility  of  the 
core  and  the  bending  rigidities  of  the  faces  in  addition 
to  the  assumption  that  faces  do  not  take  up  any  shearing 
load.  None  of  these  assumptions  have  been  made  in  the 
present  analysis.  Also,  they  have  not  considered  the  prob- 
lem of  instability.  Therefore,  excepting  the  considera- 
tions of  dynamic  effects,  the  present  analysis  can  be  con- 
sidered to  be  more  general  than  [43]. 

The  work  that  comes  closest  to  the  present  analysis 
is  the  one  by  Lisy  [47],  who  has  treated  both  the  problems 
of  bending  and  stability  of  a viscoelastic  sandwich  plate. 

The  constitutive  equations  for  the  viscoelastic  material 
considered  in  [47]  are  the  same  as  in  the  present  analyses. 
Also,  the  form  of  solution  for  w(t)  is  the  same  in  both. 
However,  as  mentioned  earlier  in  Chapter  1,  Lisy  has  made  many 
assumptions  which  have  not  been  made  in  the  present  analysis. 
Also  Lisy  has  calculated  the  critical  load  only  at  t = 00 , 
while  the  present  analysis  gives  expressions  for  the  criti- 
cal loads  for  all  t _>  0 . Additionally,  he  has  considered 


only  antisymmetric  deflection  and  buckling  while  the 
present  analysis  considers  deflection  and  buckling  for 
both  the • anti symmetri c and  symmetric  modes.  The  present 
analysis  is  therefore  much  more  general  than  the  one  by 
Li  sy . 


CHAPTER  5 


RESULTS,  DISCUSSION  AND  CONCLUSIONS 
Resul ts 

The  numerical  values  of  w(t)  and  W(t),  which  are 
the  time  dependent  antisymmetric  and  symmetric  deflec- 
tions respectively  of  a viscoelastic  sandwich  plate  con- 
sidered in  the  last  chapter,  have  been  calculated  and  re- 
ported in  this  chapter.  The  values  of  r ( a ) ( t ) anci 

N , \(t),  which  are  the  time  dependent  critical  values 
cr  ( s ) 

of  the  axial  load  for  antisymmetric  and  symmetric  modes 
of  buckling  respectively,  have  also  been  calculated. 

Since  the  present  analysis  aims  at  studying  the 
viscoelastic  behavior  of  a sandwich  plate,  each  of  the 
viscoelastic  parameters  (except  Cg ) has  been  considered 
over  a range  so  that  the  changes  in  the  response  of  the 
plate  for  the  corresponding  changes  in  these  parameters 
can  be  determined,  and  the  other  geometrical  parameters 
have  been  given  some  specific  values  which  have  been  kept 
unchanged  throughout  all  the  calculations.  All  the  quan- 
tities which  appear  in  the  calculations  and  results  are 


1 1 1 


1 1 2 


in  non-dimensional  forms.  The  following  non-dimensional 
quantities  are  defin-ed  as: 


C1 

b/a  , 

( 

c g = -/a  , 

c3  = vc2/vcl 

C4 

/b 

cl  ’ c5 

= pf2/ucl  , 

c 6 = (P-i  + P2)/^ 

C 7 

(Pi  - 

P2 ) /y c 1 5 

C8 

= nf/nc  , 

Ct  = 

t/Tc 

(where  tc 

= V'cl1 

C N 

N i / ( U 

cl  a)  ’ 

Cx  = x/a 

and  Cy  = y/b 

In  the  specific  problems  considered  in  this  chap- 
ter, the  following  non-dimensional  parameters  have  the 
values: 

v j = . 28,  v = . 3 , r = 0.067 

f c 

c-|  = 1.0  , c 2 = 0.02  , and  Cg  = 100 

-4  - 3 

In  addition,  c ^ = 0.4  x 10  , c g*  = 0.5  x 10  and  = 

_ 2 

0.1  x 10  for  the  problems  of  deflection.  Also,  all  the 

deflections  have  been  calculated  at  c = 0.5  and  c = 

x y 

0.5,  i.e.,  at  the  center  of  the  plate.  Each  of  the  param- 
eters Cg,  c ^ and  c g has  been  given  different  values  over 
a range  and  both  the  antisymmetric  and  symmetric  deflections 
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have  been  calculated  for  o <_  t < «>.  Figures  5-9  show 
these  results.  From  Figure  9,  it  is  found  that  the  sym- 
metric deflection  is  quite  insensitive  to  changes  in  the 
value  of  Cg  over  the  range  from  50  to  1000  for  the  given 
values  of  and  shown  in  the  same  figure.  Almost 
exactly  the  same  graph  as  shown  in  Figure  9 is  obtained 
if  the  value  of  C5  is  kept  fixed  at  100,  and  C4  is  al- 
lowed to  take  up  any  value  in  the  range  from  50  to  1000. 
Thus,  it  is  found  that  the  anti symmetri c deflection  at 
any  given  time  is  virtually  unaffected  by  any  changes 
in  the  values  of  C4  and  if  these  two  quantities  have 
relatively  high  values  compared  to  Cv  which  is  usually 
the  case. 

In  the  buckling  problems,  it  is  found  that  for 
each  of  the  antisymmetric  and  symmetric  modes,  the  mode 
numbers  m and  n,  equal  to  the  numbers  of  half  sine  waves 
in  the  x arid  y directions  respectively,  which  give  rise 
to  the  minimum  values  of  the  critical  axial  load,  remain 
the  same  for  0 <_  t £ °°.  Therefore,  the  values  of  m and 
n for  the  minimum  values  of  the  critical  loads  can  be 
determined  by  comparing  the  different  values  of  the  crit- 
ical loads  correspond!- ng  to  different  values  of  m and  n 
at  any  value  of  t between  zero  and  infinity.  Here,  this 
has  been  done  at  t = °°.  Once  m and  n are  determined,  the 
corresponding  critical  loads  are  plotted  against  time  for 


0 < t < “ and  shown  in  Figures  10  and  12. 

It  can  be  shown  [42]  that  both  N , , and  N , x 

are  monotone  increasing  functions  of  n;  thus,  the  minimum 

values  of  and  x occur  at  n = 1 . But  the  same 

is  not  necessarily  true  for  m.  For  antisymmetric  buckling 

it  is  found  that  for  the  particular  problem  considered 

here,  the  minimum  value  of  critical  axial  load  occurs 

either  between  m = 1 and  m = 2,  or  at  m = 1;  and  in  the 

former  case,  the  values  at  m = 1 and  m = 2 differ  very 

slightly  from  each  other  (the  value  at  m = 1 being  the 

smaller  one).  Therefore  N , » has  been  calculated  for 

c r \ a ) 

m = 1 and  for  0 < t < “. 

For  symmetric  buckling,  the  value  of  m correspond- 
ing to  the  minimum  critical  load  is  very  high.  This 

1 s shown  in  Figure  11. 

In  the  calculations,  the  following  non-dimensional 
parameters  denoting  the  critical  loads  have  been  used: 


Cna  * Nc(a)/Kla> 


and 


C = N , \ / ( u t a ) 
ns  cr(s)VMcl  ' 


Figures  5 through  12  show  the  results. 
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Figure  5.  Anti symmetri c deflection  for  different  values  of  vs  time 
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Figure  6.  Symmetric  deflection  for  different  values  of  c3  vs  time 


Figure  7.  Anti symmetri c deflection  for  different  values  of  c„  vs  time 
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Figure  8.  Antisymmetric  deflection  for  different 
values  of  c ^ vs  time 


Figure  9.  Symmetric  deflection  for  50  < c5  < 1000 
vs  time 
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Figure  10.  Critical  load  for  antisymmetric  buckling  vs  time 
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Figure  11.  Determination  of  mode  number  for  symmetric 
b u c k 1 i n g 
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Figure  12.  Critical  load  for  symmetric  buckling  vs  time 
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Discussion 

A viscoelastic  material  represented  by  a three- 
parameter  model  shown  in  Figure  4 in  a quasi-static 
problem  can  be  interpreted  as  an  elastic  material  whose 
elastic  modulus  is  changing  with  time,  from  that  of  a 
single  spring  of  stiffness  equal  to  y c2  (if  the  material 
considered  is  that  of  the  core)  at  t = 0,  to  that  of  a 
combination  of  two  springs  of  stiffnesses  yc2  and  y ^ 
in  series,  at  t = At  any  other  time  t where  0 < t<  °° , 

the  equivalent  elastic  modulus  is  that  of  a combination 
of  springs  of  stiffnesses  uc2  and  K(t)ycl  in  series,  the 
value  of  the  factor  K(t)  being  determined  by  q , and 
K(0)  = oo,  k(°°)  = 1.  Similar  interpretation  holds  for  the 
material  of  the  faces. 

Thus  it  becomes  apparent  that  the  response  of  the 
viscoelastic  sandwich  plate  to  a given  set  of  loading  can 
be  obtained  directly  from  that  of  an  elastic  plate  in  two 
special  cases  viz.  at  t = 0 and  at  t = These  results 

can  be  used  as  a check  on  the  results  obtained  from  the 

analysis  of  viscoelastic  sandwich  plates.  This  is  done 

% 

in  the  following: 

Symmetric  Deflection 


For  the  symmetric  deflection,  it  is  readily  seen 


from  Equations  (139)  and  (140)  of  the  viscoelastic  analy- 
sis that 


“mo'*  ■ 0)  ' H45 


and 


W (t  = * ) = 
mn 


H 4 5 H 4 4 


‘47 


(157) 


From  (138),  (112)  and  (120)  through  (123),  it  is 
found  that  is  equal  to  the  value  of  W as  given  in 

Equation  (111)  of  the  elastic  analysis  if  Gc  and  are 
replaced  by  yc2  and  yf2  respectively.  Also  (H45H44)/H47 
is  equal  to  the  value  of  W as  given  in  Equation  (111) 
of  the  elastic  analysis  i f Gc  and  Gf  are  replaced  by  Kce 

and  K_pe  respectively,  where 


K = yc2  ycl 
ce  y , y 

c2  cl 


stiffness  of  the  combination  of  two 
springs  of  stiffnesses  y ? and  y , in 
series  c c 


and 


„ _ yf2  yfl 

f e y , y 

f 2 f 1 


stiffness  of  the  combination  two  springs 
of  stiffnesses  y ^ 2 and  y ^ -|  in  series 


Anti symmetri c Deflection 


As  for  the  antisymmetric  deflection,  it  is  found 
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from  Equation  (141)  that 


“mn(t  V“>  = H30B1 


(158) 


and  from  (125)  and  (120)  through  (123)  it  is  seen  that 

H 3 0 B l eclua^  to  the  value  of  as  given  in  Equation 

(110)  of  the  elastic  analysis  if  Gc  and  Gf  are  replaced 

by  Kce  and  K p respectively. 

To  find  the  value  of  w (t  = 0),  a direct  way 

mn 

would  be  to  put  t = 0 in  one  of  the  Equations  (129) 
through  (136).  But  this  involves  very  long  algebraic 
expressions  which  are  difficult  to  work  with.  To  get 
around  this  difficulty,  the  following  procedure  can  be 
f o 1 1 owed : 

If  f(s)  = Laplace  transform  of  F(t), 

then 

Lt  s f ( s ) = F ( 0 ) (159) 

S-J-oo 


Using  (159)  and  (124)  the  following  is  obtained: 


w ( t = 0 ) = 
mn ' ' 


Lt  sw. 

S->oo 


mn 


(160) 
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It  is  now  readily  seen  from  (125)  and  (110)  that 
H30  s egual  t0  t*ie  va^ue  wnm  as  given  in  Equation  (110) 
of  the  elastic  analysis  if  G and  Gf  are  replaced  by  y 
and  respectively. 

An  independent  way  of  calculating  the  critical 
loads  Ncr(a)  and  Ncr(s)  at  t = 0 and  at  t = <*>  is  to  con- 
sider the  expressions  wmn(t  = 0),  wmn(t  = ®) , Wmn(t  = 0) 

and  Wmn  ^ t = 00 ) as  given  in  ( 1 60 ) , (1  58)  and  (157).  This 
can  serve  as  a check  on  the  results  obtained  in  the  last 
chapter  for  wmn(t)  and  W (t)  where  0 < t < ».  This  is 
done  as  foil ows : 

Antisymmetric  Buckling 

From  (160)  and  (125),  it  is  found  that 

H, 

"mnO  = °>  = H30  = ~~ 


The  only  way  in  which  wmn ( t = 0)  can  become  infinity  for 
finite  transverse  loads  is  when  H25  goes  to  zero,  and  the 
condition  H25  = 0 gives,  from  (12,5), 


Nn.(t  = 0) 


Ncr(a) ^ = °) 


Pc6  + H18Pc6Pf6  + H19Pf6 
H20Pc6  + H21 Pf6 


(161) 


Also, 
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[using  (159)] 


H 


56 


(162) 


From  (161)  and  (162),  it  is  found  that  the  same  expression 
for  Ncr(a)(t  = 0)  is  obtained  in  two  different  ways. 


(t  = °°)  is  calculated  first  from  w (t  = °°)  and  then  from 

mn 

the  expression  for  Ncr(a\(t),  the  same  value  is  obtained 
as  given  in  (158). 

Symmetri c Buckling 

Each  of  the  values  of  N , N(t  = 0)  and  N (t  = °°) 

cr  ( s ) cr 

can  be  obtained  in  two  different  ways  and  can  be  shown  to 
be  the  same  following  a procedure  similar  to  that  for  anti- 
symmetri c buckling. 

One  important  point  is  to  be  noted  here.  The  value 

of  the  critical  load  N,\  can  be  obtained  from  the  ex- 

c r ^ a ) 

pression  for  w (t)  as  shown  in  this  section  only  at  t = 0 
and  t = °°.  At  any  other  value  of  t where  0 < t < <*>,  it  can 
be  shown  that  no  definite  condition  can  be  obtained  for 
which  w (t)  becomes  infinity.  So,  to  obtain  the  value  of 
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Ncr(a)(t)  for  any  time  t where  0 < t < °°,  one  must  cal- 
culate it  from  Equations  (129)  through  (136). 

The  same  is  true  of  N , >t. 

c r ( s ) 

Applicability  of  the  Present  Analysis 

The  present  analysis  enables  one  to  predict  the 
deflection,  antisymmetric  or  symmetric,  at  any  time  after 
the  application  of  loads,  of  a viscoelastic  sandwich 
plate.  The  more  practically  useful  aspect  is  the  buckling 
phenomenon.  The  present  analysis  also  predicts  the  buck- 
ling phenomenon  for  0 £ t £ °°.  Figure  13  shows  a graph 
(like  the  ones  obtained  in  the  present  analysis)  of  a crit- 
ical axial  load  for  antisymmetric  or  symmetric  mode  of 
buckling,  vs  time  t for  0 £ t £ °°  in  a given  problem. 

The  curve  ABC  shows  the  critical  value  of  the  axial  load 
N (t).  The  region  above  the  curve  ABC  is  unstable,  and 
the  region  below  ABC  is  stable. 

Suppose  at  t = 0,  the  plate  is  subjected  to  a 
step  load  which  is  axial  (in  the  x direction)  and  compres- 
sive. If  the  magnitude  of  this  load  is  smaller  than 
N (t  = °°)  represented  by  the  line  D'B',  the  whole  line 
D'B'  falls  into  the  stable  region,  and  therefore  the  plate 
will  not  buckle  for  any  time  t,  0 £ t £ If  the  magni- 
tude of  the  load  is  greater  than  fJ  ( t = 0),  the  line  D"B" 
representing  this  case  is  completely  in  the  unstable  region, 
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Figure  13.  Applicability  of  the  analysis  of  buckling 
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and  the  plate  buckles  immediately  after  the  application 
of  the  load.  Finally,  if  the  magnitude  of  the  load  is 
greater  than  N£r(t  = «)  but  less  than  N (t  = 0),  this 
is  shown  by  the  line  DB  which  is  in  the  stable  region 
for  0 _<  t < t1,  and  in  the  unstable  region  for  t'  < t < °°. 
Therefore  buckling  occurs  at  t - t'. 

Therefore  an  important  prediction  can  be  made  that 
the  plate  will  buckle  at  t = t 1 after  application  of  the 

load.  Conversely,  if  the  plate  is  required  to  stand  a 

load  for  only  a given  length  of  time,  say  t 1 , the  corre- 
sponding maximum  load  N = N (t  = t')  can  be  calculated 

in  a reverse  way. 

Conclusions 

The  analysis  presented  in  this  dissertation  con- 
sists of  two  main  parts  connected  by  a comparatively  short 
but  important  link  in  the  middle.  The  first  of  the  main 
parts  is  a nonlinear  theory  of  elastic  sandwich  plates 
derived  in  a unified  approach  from  the  variational  prin- 
ciple of  Hellinger  and  Reissner  in  three-dimensional  elas- 
ticity. It  can  be  considered  to  be  an  extension  of  the 
earlier  wo r ks  o f Ebc i og 1 u and  Habip.  Ebcioglu  considered 
a more  general  problem  (of  a sandwich  plate),  whereas  Habip 
considered  a single  layer  plate.  Habip  used  a more 

general  procedure  (from  the  static  and  quasi-static  elasticity 
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point  of  view)  viz.,  the  Modified  Hellinger-Reissner 
principle,  than  Hamilton's  principle  which  was  used  by 
Ebcioglu.  The  present  analysis  has  used  Habip's  pro- 
cedure for  the  problem  considered  by  Ebcioglu.  Kim's 
linear  analysis  for  the  problem  of  a sandwich  plate  sub- 
jected to  combined  axial  and  transverse  loads  has  been 
derived  as  a special  case  of  the  present  analysis. 

The  second  of  the  two  main  parts  of  this  disser- 
tation is  the  solution  of  the  problem  of  a viscoelastic 
sandwich  plate  subjected  to  combined  axial  and  transverse 
loads.  This  part  is  an  extension  and  development  of  ear- 
lier analyses  of  viscoelastic  sandwich  plates  by  other 
authors.  The  present  analysis  has  used  a different  pro- 
cedure (the  principle  of  correspondence)  than  those  used 
by  others,  and  has  introduced  a new  idea  of  time  dependent 
critical  loads. 

The  main  two  parts  of  the  dissertation  described 
above  are  connected  by  the  well-known  correspondence 
principle  in  linear  viscoelasticity,  which  has  been  de- 
scribed as  the  "short  but  important  link"  in  the  beginning 
of  this  section.  The  viscoelastic  analysis  has  been  de- 
rived from  the  elastic  analysis  using  the  correspondence 
principle. 

Further  extensions  of  the  present  work  can  be 
carried  out  on  many  fronts  in  both  the  fields  of  elasticity 
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and  viscoelasticity.  For  elastic  sandwich  plates,  different 
and  more  general  types  of  loading  and  boundary  conditions 
can  be  considered.  For  viscoelastic  sandwich  plates,  a 
f our-parameter  model  or  a still  better  model  can  be  used. 
Also,  nonlinear  problems  of  viscoelastic  plates  can  be 
considered  using  the  governing  equations  of  nonlinear  vis- 
coelasticity. Finally,  on  the  numerical  side,  there  are 
and  will  always  be  plenty  of  useful  problems  on  sandwich 
plates. 
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